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ABSTRACT 

Eleven bright gamma-ray bursts (GRBs) detected by BATSE have also been seen at 
much higher energies by EGRET, six at energies above 10 MeV. Most distinctive among 
these is GRB940217, which includes long duration, hard gamma-ray emission and the 
most energetic GRB photon detection to date, around 18 GeV. Such observations imply 
that these bursts are optically thin to photon-photon pair production at all observed 
energies, for target photons both internal and external to the source. For bursts more 
than about 30pc away, internal transparency can be achieved only if the source is moving 
with a relativistic bulk Lorentz factor F ^ 1, or if the radiation is highly beamed. Early 
calculations of 77 e'^e~ considerations for GRBs were limited to cases of a beam with 
opening half-angle Gb ~ 1/F, or expansions of infinitely thin spherical shells. This paper 
presents our extension of pair production optical depth calculations in relativistically 
expanding sources to more general geometries, including shells of finite thickness and 
arbitrary opening angle. The problem is reduced analytically to a single integral in the 
special, but quite broadly applicable case, of observing photons only along the axis of the 
expansion. We find that the minimum bulk Lorentz factor for the EGRET sources to be 
optically thin, i.e. display no spectral attenuation, is only moderately dependent on the 
shell thickness and virtually independent of its opening solid angle if 0b ^ 1/F. This 
insensitivity to 0b relieves the commonly-perceived number problem for non-repeating 
sources at cosmological distances, i.e. it is not necessary to invoke small Gb to effect 
photon escape. The values of F obtained, typically of the order of 10 for halo bursts 
and ^ 100 for sources of cosmological origin, depend somewhat on the choice of GRB 
timescale used to determine the expansion size. Our new limits on required velocity for 
given source geometries will aid in placing realistic constraints on GRB source models. 
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1. INTRODUCTION 

Gamma-ray bursts (GRBs) are the brightest sources in the gamma-ray sky and they may also 
be among the most distant sources in the Universe. The discovery by the BATSE detector on the 
Compton Gamma-Ray Observatory (CGRO) that the spatial distribution of GRBs is isotropic and 
inhomogeneous (Meegan et al. 1992, 1996) suggests that the sources arc cither in an extended 
galactic halo or at cosmological distances. The level of isotropy of the GRB spatial distribution 
limits halo models to core radii of around 50-80 kpc (Hakkila et al. 1995); tighter constraints 
are expected for more recent data accumulations, and Briggs et al. (1996) suggest that a galactic 
halo shell distribution must be at least 120 kpc distant. The observed average fluxes of GRBs at 
Earth therefore imply high luminosities for isotropically emitting sources: L ~ 10^2-43 gj.g g-i 
a distance of d = 100 kpc and L ~ 10^''"^^ ergs~^ at a distance of d = 1 Gpc. In addition, rapid 
time variability (At ~ several ms) is observed in GRB light curves, whose structural diversity is 
illustrated, for example in the BATSE IB catalogue in Fishman et al. (1994). This variability 
implies a compact source size, which in combination with the high luminosities yields photon 
densities that are high enough to make galactic halo or cosmological GRBs optically thick to 
photon-photon pair production by many orders of magnitude. One would then expect attenuation 
of the observed spectrum (perhaps as a quasi-exponential cutoff, trough or shelf: examples are 
depicted in Baring and Harding 1997) around the pair production threshold of 1 MeV if the GRB 
sources are more distant than a few kpc and have quasi-isotropic radiation fields (Schmidt 1978, 
Epstein 1985). 

Yet GRB spectra arc observed to extend well beyond 1 MeV and into the GeV range. The 
GRS detector on the SMM satellite first measured emission in GRB spectra significantly above 1 
MeV, often extending up to 10 MeV (e.g. Nolan et al. 1983), and in one case up to 80 MeV (Share 
et al. 1986). BATSE routinely observes GRB spectra extending up to and above 1 MeV. While 
most bursts exhibit spectral steepening at a variety of energies between 50 keV and few hundred 
keV (Band et al. 1993; see also Schaefer et al. 1994 for the BATSE IB spectroscopy catalogue), 
a number of bursts display spectral breaks between 500 keV and about 2 MeV (Schaefer et al. 
1992), but no cutoffs. The EGRET instrument, also on CGRO, has detected emission above 50 
MeV from four of the brighter GRBs triggered by BATSE, a fifth up to 30 MeV and another three 
up to a few MeV; all are consistent with power- law spectra extending to as high as 1.2 GeV, in the 
case of GRB930131 (Sommer et al. 1994), and 3.4 GeV for GRB940217 (Hurley et al. 1994). The 
GRB940217 source is best known for exhibiting delayed or prolonged high energy emission, detected 
80-100 minutes (i.e. more than one full earth orbit of CGRO) after the initial trigger, including a 
photon of energy 18 GeV (Hurley et al. 1994) that is not markedly inconsistent with extrapolation 
of the power-law continuum. In fact, some evidence for delayed high energy emission pre-dated 
GRB940217, with the observation (Dingus et al. 1994) of a single 10 GeV photon that could have 
been associated with GRB910503. It is clear that, in contrast to soft gamma-ray repeaters (SGRs), 
which are now believed to probably be a separate class of galactic sources (although the classical 
GRB behaviour reported by Fenimore, Klebesadel and Laros 1996 of the 5th March 1979 outburst 
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at early times can support proponents of GRB-SGR associations), there have been no attenuation- 
type turnovers or cutoffs observed in a GRB spectrum. High energy gamma-ray emission therefore 
may be common in bursts, and the EGRET detection rate is consistent (Dingus 1995, though this 
inference is subject to poor statistics) with ah bursts emitting above about 30 MeV. Observed GRB 
spectra are therefore in direct conflict with predicted pair production cutoffs for isotropic emission. 

An obvious solution (e.g. Krolik and Pier 1991, Fenimore, Epstein and Ho 1992) is to allow 
some anisotropy of the emission, so that the interaction angles Ou of the photons are restricted. 
The threshold for pair production, St = 2/[£j(l — cos^)] where et and £j are the energies of a 
test photon and an interacting photon in units of mgC^, could therefore be increased above the 
maximum observed energy. Beaming of the radiation can be achieved through relativistic motion: 
the radiation from a source that is isotropically emitting in the comoving frame will be beamed 
mostly within an angle of the order of 1/T in the observer's frame, where T is the bulk Lorentz 
factor. For the case of a small emitting blob moving relativistically, the pair production optical 
depth r-y-y is reduced by a factor r^^^^^") below the optical depth for isotropic radiation, where a 
is the photon spectral index (Krolik &; Pier 1991, Baring 1993). The minimum bulk Lorentz factors 
required to make t^j < 1 in the bright "superbowl" burst (GRB930131) detected by EGRET 
(Sommer et al. 1994) up to an energy of ~ 1 GeV are F ^ 10'^ at a distance of 1 Gpc and F ^ 10 
at 30 kpc (Harding 1994, Harding & Baring 1994). In this case of relativistic beaming within angle 
1/F, the required luminosity L at the source is smaller because the observed flux, ~ r'^L/Aircf, is 
enhanced by a solid angle factor F^ (Krolik & Pier 1991). However, the number of sources must be 
a factor F^ higher in order to account for the observed number of GRBs. In the case of cosmological 
GRBs, this factor could be as high as 10^ for the above limits on F. This is unacceptably large 
for many of the proposed models, including neutron star-neutron star or neutron star-black hole 
mergers (Paczyhski 1986; Eichler et al. 1989; Narayan, Piran and Shemi 1991; Meszaros &: Rees 
1992), failed Type lb supernovae (Woosley 1993) and rapid spin-down of high-field millisecond 
pulsars (Usov 1992), and hence defines the so-called "number problem" for beamed cosmological 
bursts. 

Source geometries with beaming angles larger than 1/F could ease this problem if the high 
energy photons were able to escape. In fact, the radiation from GRB sources in the Galactic 
halo or at cosmological distances is expected to involve a wind or fireball expanding relativistically 
(Paczynski 1986, Goodman 1986, Piran and Shemi 1993). Fenimore, Epstein & Ho (1993) have 
shown that a relativistically-expanding, thin spherical shell will allow escape of high energy gamma- 
rays, because a test photon on the surface of the shell will not be able to interact with all other 
emitted photons due to causality limits. This arises as a consequence of the transient nature of the 
emission, since then only photons emitted within a "look-back" surface around the test photon will 
interact to contribute to the pair production optical depth. The Fenimore, Epstein &; Ho (1993) 
calculation was limited to the case of an infinitely thin shell. 

In this study, we have extended the calculation of the pair production optical depth in GRB 
sources to the full range of source geometries: opening angles from 1/F to a spherical expansion. 
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and shells of arbitrary thickness. The optical depth for test photons emitted within the expanding 
shell will be limited to interaction with other photons within a "look-back" volume. We present an- 
alytic development and simplification of the pair production optical depth, make detailed numerical 
calculations, and derive analytic expressions in various limits. The intent is to provide a model- 
independent evaluation of the pair production opacity of a relativistically expanding, transient 
gamma-ray source whose emitted spectrum extends above observed energies. Using our results, we 
derive estimates for the minimum bulk Lorentz factors required for source transparency in those 
GRBs detected by EGRET at high energies. These limits are largely insensitive to the source 
opening angle provided that it exceeds 1 /T . This reflects the strong impact of causality in deter- 
mining the optical depth, and clearly renders the number problem for cosmological source models a 
non-issue: the total negation of this number problem for a wide range of expansion geometries is a 
principal conclusion of this paper. A detailed description of the source geometry and the derivation 
of an analytic form for the pair production optical depth (and associated limiting cases) for infinite 
power-law source spectra are presented in Section 2; there the general quintuple integral expres- 
sion for Eq. (pO[)] is expediently reduced to a comparatively simple single integration in 



Eq. (37), a principal result of this research, rendering our developments quite amenable to various 
observational and theoretical applications. Section 3 is devoted to the application of these results 
to EGRET bursts and the estimation of minimum bulk Lorentz factors in these sources, including a 
discussion of the behaviour of our results in relevant parameter spaces and various issues pertaining 
to our calculations. Readers more interested in the applications and implications of our calculations 
than in the detailed derivations presented in Section 2 should note that Eq. (37) is the final form 
for the pair production optical depth, which should be used in conjunction with the normalization 
specified via the flux in Eq. (|l8|) . 



2. THE PAIR PRODUCTION OPTICAL DEPTH 

The generic picture of a gamma-ray burst that is considered here assumes the photon source 
(i.e. region of emission) to be expanding with constant and homogeneous bulk Lorentz factor 
r , with opening angle 20b about some axis and thickness AR (constant throughout) in the 
observer's frame, and with the initial condition (at time t = 0) that the source's inner radius 
is -Ro ■ The expansion therefore traces a conical volume and can assume a variety of geometries 
such as solid cones, solid spheres or spherical shells, at any given time, depending on the values of 
the input parameters Gb and AR/Rq . The constancy of T in time is a convenient assumption 
that is not strictly valid during early epochs of the expansion if AR/Rq is not much less than 
unity: when AR/Rq ^ 1 , the expansion initially resembles a quasi- isotropic and almost stationary 
radiation gas in the observer's frame. The adiabatic redistribution of momenta that naturally 
occurs in expanding (and "inert") photon gases is therefore neglected. We also opt to ignore the 
consideration of possible dynamic acceleration or deceleration of the underlying plasma, since such 
dynamics are quite model-dependent. 
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Suppose that a test (i.e. potentially observable) photon, with energy etmeC^ , is emitted at 
time t = from the inner radius Rq of the shell and moves through the source to eventually 
escape and reach the observer. We opt for test photons originating at the back of the expansion 
throughout this paper; starting these photons closer to the outer surface will reduce the optical 
depth they encounter by a factor of order unity, so that the results we obtain will be qualitatively 
representative of general initial positions for test photons within the source. If the angle cosine 
between the test photon's momentum and position vectors is fit , then the radial distance rj of the 
test photon from the center of the expansion, at any time t , is 

rt = Ro + iJ-tct . (1) 

The overall geometry of this source expansion is depicted in Fig. ^ and is discussed in more 
detail below. Radial propagation of the test photon along the axis of the expansion corresponds 
to /if = 1 . An important assumption about the expansion that is made in this paper is that 
no photons are present prior to time t = . This "switch-on" stipulation restricts the photon 
population that can causally interact with the test photon at early times, and indeed mimics burst 
temporal behaviour; it is anticipated that the details of the switch-on will have only a quantitative 
rather than a qualitative influence on the results presented. The objective of this section is to derive 
an analytic expression for the pair production optical depth for this expanding source geometry. 
Consideration of the influence of the plasma that is present in the emission region on the photons 
it generates will be omitted from this analysis. 

The optical depth for two-photon pair production 77 — > e^e~ can be obtained from well- 
known expressions for the reaction rate R^^ for interactions of photons in a single population (e.g. 
see Eq. (27) of Weaver 1976, or Eq. (7) of Stepney and Guilbert 1983). In the case where one 
photon is a test photon of dimensionless energy et , the optical depth, differential in the distance 
rt that the test photon travels, is (e.g. see Eq. (7) of Gould and Schreder 1967) 

'^^dri^*'' ^ ^/ '^77(x) (1 -^ti)^(ei; /"ti; ^t)c^ei^^Mti • (2) 

Here subscripts i denote quantities of the photon that interacts with the test photon, x = 
^/et£i{l — Uti)/'^ is the center-of-momentum (CM) frame energy scaled by m-eC^ , and 

^77(X) = ^[[2x^ + 2x^-1) \ogJyx+^X^-l]-x{l+X^)^X^-l] ■ (3) 

is the Lorentz-invariant pair production cross-section (e.g. see Eq. (13-40) of Jauch and Rohrlich, 
1980), where r^. = e^/nieC^ is the classical electron radius. Hereafter, all photon energies will 
be assumed to be dimensionless, being scaled by mgC^ . Also, ^xti = cos Ou is the angle between 
the momentum vectors of the test and interacting photons, and n(ei, ^lu', i~t) is the source photon 
density distribution at the position of the test photon. The factor of 1/2 in Eq. (^) is the standard 
correction for double-counting in interactions of identical particles; it is omitted in the calculations 
of Epstein (1985; see his Eq. 2.8) and Zdziarski (1984), who treat the test photons as a separate 
population from the interacting photons. 
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It is instructive to identify the typical energy Ei of photons that interact with test photons at 
Et , specificahy for an expansion of bulk Lorentz factor T that spawns power-law photon spectra, 
the conditions pertaining to the analysis of this paper. For such test photon energies, the minimum 
possible energy of the interacting photons, defined by the pair production threshold, is ~ r/e^ in 
the rest frame of the expansion, and of the order of V"^ /et in the stationary observer's reference 
frame. The pair production cross-section in Eq. (|3|) peaks not far above threshold, and since the 
optical depth is a convolution of this cross-section and the spectrum, which is a strongly decreasing 
function of energy, it is clear that the typical energy of an interacting photon is usually never far 
above threshold, i.e. around V"^ /et in the observer's frame. This result holds regardless of the 
source photon density provided phase space near threshold is accessible, which always is the case 
for infinite power-law spectra. 



2.1. Source Geometry 

Before developing the calculation of the pair production optical depth it is both necessary and 
elucidating to elaborate the details of the expansion geometry and define useful spatial variables. 
The general picture of the expansion at any instant is given in Fig. [l|a; the definitions of the test and 
interacting photons' spatial and angular variables are depicted in Fig. ^3 and are now enunciated. 
The radius of the test photon at any time is given by Eq. (||) and the angle between the radius 
vector (OT in Fig. |^b) to the test photon and the cone axis (OZ in Fig. ^a), which bisects 20b , is 
defined as Qt ( ^ ©b ) • If the expansion has an inner radius R at time t and has a thickness Ai? 
that is constant in time, then the test photon remains within the expanding volume only when 

R = Rq + (ict < n < Ro + AR + Pct , (4) 

where /3 = \/l — 1/T'-^ is the bulk velocity of the expansion (in units of c). In general, the 
angle between the test photon's position and momentum vectors is 9t = arccos /it ; however, unless 
otherwise stated, in subsection 2.1 and subsequent portions of the paper, 9t will be assumed to be 
zero so that rt = Rq + ct ; this will be a specialization to the most salient case of radial propagation 
of test photons. Eq. (^) leads to the determination of the time te the test photon takes to escape 
the expanding plasma: 

= c(l - /?) ■ 

Of course, calculation of the pair production optical depth will involve an integration over times 
< t < oo , including when the test photon has escaped the expanding plasma. 

The test photon interacts with photons at positions within some causally-connected look-back 
volume. Detailed considerations of such look-back regions for relativistic expansions are presented 
in Rees (1966), Fenimore, Epstein and Ho (1993). Suppose that a typical interacting photon is 
located at a radius rj with the angle between its position vector (01) and the expansion axis (OZ) 
being Oj (see Fig. ||b). Such an interacting photon was emitted at time U {0 < ti < t) and at 
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a distance ru = c{t — ti) from the test photon. Let Qu be the angle between the radius vectors 
(OT and 01 in Fig. |l]b) of the test and interacting photons. Further let the angle between the 
momentum vectors of the test and interacting photons be 9^ . If the line (TI in Fig. |^b) between 
the positions of the test and interacting photons makes an angle 6i with the radius vector of the 
interacting photon, then simple geometrical analysis gives 

= + - '^nri cos Qti , 
rf = + rfi - 2rtrti cos{0i - @ti) , (6) 
rt = rf + r^j + 2rirti cos 0i . 

These relationships will be used to develop the integrations in the expression for the optical depth 
that is derived in the following subsection. In general, if (p is the angle between the planes defined 
by the test photon momentum and position vectors and the momentum vectors of the test and 
interacting photons, then spherical trigonometry yields 

cos(^j — Qti) = cos 6t cos 6ti — sin 9t sin Ou cos (p . (7) 

However, when specializing to the case of radial propagation of test photons, the photon momenta 
lie in the OTI plane so that (p = ir and 9i = Ou + Qu ; this simplification will be used in subsequent 
sections. 

The geometry of the source defined above, restricts the values of the variables ru and Ou that 
prescribe the position of the interacting photons. These restrictions arise because the interacting 
photons can only be emitted from the portion of the region that the expanding plasma occupied 
at the time of emission that is causally connected to the test photon. In the radial direction, the 
volume that the emission region occupies at time t = ti is specified simply by 

R^ < Vi < R^ + AR , R, = Ro + pcti , (8) 

where the relation ru = c{t — ti) implies that i?j = (1 — P)Ro + P{rt — ru) . With the aid of Eq. (|^, 
this radial constraint becomes 

rf + rl - {Ri + AR)^ . m n ^ ^ r,^ + rj - Rf 

which, for the case of radially-propagating test photons (i.e. when 9i — Qu = Ou), is a compact 
representation of the limits to the nu integration in Eq. (|2|) . The values of ru that are achievable 
are further constrained by the causality condition 

< ru < ct = ^^^^ . (10) 

In fact, this restriction automatically guarantees that {rf + rf^ — Rf) / {2rtru) > 1 and therefore that 
the right hand inequality of Eq. (|9|) is always satisfied. Physically this occurs because interacting 
photons emitted at the rear of the expanding volume can never catch the test photon. 
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The constraints that the emission volume places on angles are simply enunciated. The require- 
ment that the interacting photon be within the expanding shell imposes no restriction on the radial 
variables, but does constrain 0j according to 

< < Gb . (11) 

A similar condition limits the values of at t = . Eq. (|ll]) restricts the allowable azimuthal 
angles of the interacting photon for off-axis propagation of test photons; discussion of this restriction 
is deferred to section 2.2.2 below, specifically Eq. (p5|). In addition, the range of Qu is clearly 
bounded by the expansion geometry. Inspection of Fig. ||b reveals that the maximum possible 
value of Qti is Gb -|- Gt ; with the aid of Eq. (0) this becomes the global Qu constraint 

cosGti = . '^-P^^t^ > ^ = cos(min{7r, Gb + Gj) , (12) 

y 1 - 2pfiu + 

where p = ru/rt is a scaling of ru that proves convenient in the algebraic manipulations of this 
paper (see immediately below). This can be inverted to find the ranges of acceptable values of 
pti = cos Oti , as is outlined below in Eq. (^). This concludes the presentation of the general forms 
for the constraints the source geometry places on the spatial variables defined; specific developments 
in subsequent sections are made according to algebraic need. 

It shall prove convenient to define three dimensionless variables that will facilitate the algebraic 
developments of the optical depth that are performed in this paper: 

= 7? ' ^ = ' ^ = T^P ■ (13) 

The first two of these are scaling transformations that define the coupling of length scales in the 
expansion; they are used to reduce the number of integrations in the optical depth over spatial 
variables by one [see Eq. (p^)]. Note that the causality condition in Eq. (|^) yields < p < I for 
Ht = 1 . The definition of s is effectively an alternative to the angular variable pu that proves 
convenient in reducing the integration of interacting photon angles analytically (see Section 2.3). 
These three variables will be referred to extensively in subsequent equations. 

Hereafter, this paper addresses the special case of radial propagation of test photons along 
the axis of the expansion, so that pt = cos 9t = 1 and G^ = . While this choice is motivated 
by the simplifications it introduces to the analysis, it is concordant with the goal of obtaining 
representative estimates of source bulk Lorentz factors that are consistent with GRB observations. 
Relativistically expanding sources contribute most of their observable emission along the direction 
of motion, corresponding to radial and on-axis propagation of the test photons. Off-axis (i.e. non- 
radial) emission will mostly be outside the peak of beamed radiation, and therefore form only a 
minor part of the observable flux of gamma-rays. Hence we expect that photons produced somewhat 
off-axis will contribute minimally to the observed flux, and therefore be largely irrelevant to the 
determination of minimum bulk Lorentz factors. A brief discussion of this specialization, in the 
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light of the results obtained, is presented in Section 4. Note that while "limb" photons move on 
average at large angles to the beamed photon population in their locale, the phase space that 
connects them causally to the remaining photon population is small; it is not clear whether or not 
limb photons will have enhanced pair production optical depths relative to line-of-sight radiation. 

2.2. Analytic Reduction of The Optical Depth 

The differential optical depth in Eq. (Q) can be developed once the photon distribution function 
n(ej, fiti] rt) is known, and eventually an integration over rt will be performed to obtain the total 
optical depth r-y^(et) . The form that n(ej, fiu] rt) takes depends on basic assumptions about the 
expanding photon gas. In this paper, the rate of photon emission is taken to be constant in time and 
space after time t = (following, for example, Fenimore et al. 1993, but in contrast to the uniform 
density choice made by Harding and Baring 1994), but is zero for t < 0, and is isotropic only in 
the comoving frame of the expansion. This simplifying assumption is made for its convenience; it is 
unrealistic since it may be acausal for some initial conditions (e.g. the instantaneous "switch-on" 
over a finite volume). Observed temporal behaviour in individual sources is somewhat chaotic (e.g. 
see the BATSE IB catalogue, Fishman et al. 1994), so the assumption of constant emissivity is not 
truly accurate, depending on the timescale of specific consideration. However, in conjunction with 
the source geometry prescribed here, a constant emissivity may be able to approximately produce 
global properties of bursts, such as longer average decay timescales than rise times (Nemiroff et al., 
1994; Norris, et al. 1994; Mitrofanov 1995). Our assumption of a temporally-invariant emissivity 
after switch-on is an appropriate approximation for the objectives of this paper, since only estimates 
of the bulk relativistic motion in gamma-ray bursts are at present obtainable as the origin of bursts 
is still uncertain. Note that the rate of emission might be expected to decline in time due to 
expansion effects such as adiabatic cooling and a decrease in plasma density. 

An immediate consequence of this approximation is that the photon distribution in the co- 
moving frame of the expansion is anisotropic even for isotropic photon emission, due to radiative 
transfer effects in finite source volumes. This assumption is clearly different from the premise of 
Gould and Schreder (1967), from which the bulk motion analyses of Krolik and Pier (1991) and 
Baring (1993) are derived, who all effectively assumed that the photon distribution is everywhere 
and at all times isotropic (Krolik and Pier 1991 actually invoked the equivalent assumption of 
isotropy of the photon intensity in the comoving frame). Isotropy of the photon distribution is an 
even more elementary approximation: it is perhaps less realistic than the assumption of isotropic 
emission rate that is made here, given that isotropic radiation fields are usually best generated 
in optically (Thomson) thick media, contrary to the basis of these pair-production calculations. 
Isotropic emissivities require only that the supporting particle population is isotropic, and that 
there is no other preferred direction in the emission region, such as that imposed by the presence 
of a magnetic field. 
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2.2.1. Photon density and the observed flux 

For a constant rate of photon emission, the photon distribution at any point is simply obtained 
by adapting standard radiative transfer results (e.g. see Rybicki and Lightman, 1979) to the 
consideration of photon densities. The conservation of photon numbers and volume elements along 
light rays in optically thin environments automatically yields conservation of the number density, 
from which the transfer equation yields 

n{ei, fitf, rt) = ^ J h{£i, ^f, ri) dm . (14) 

Here ru traces out the photon path through the look-back volume and n(ei, /Uj; rj)/27r is the rate 
of photon emission in the observer's frame, per unit solid angle, at the position of the interacting 
photon (the point labelled I in Fig. |l]b). The arguments of the emission rate in the integrand are 
implicitly functions of the spatial variables relating to the test photon position, i.e. Hi = cos 9i = 
Aii(n, tti, tJ-ti) and = rj(rf, tu, Hu) are defined by Eq. (0), and determined by the geometry in 
Fig. Here = cos Ou . Throughout the following analysis, it is assumed that 0^ = . The 
angular dependence of the distribution and emission rate has been retained because of the highly 
anisotropic conditions encountered in this calculation. The azimuthal angle (f)ti is defined to be the 
angle between the OZT and OZI planes in Fig. |l|b. Since the emission rate is independent of position 
within the expanding volume, no azimuthal dependence appears in the arguments of n(ej, /i^; rj) . 
Note that Harding and Baring (1994) combine Eqs. (|2|) and ( p^ ) into a single expression for the 
optical depth in their Eq. 2; their result is mildly erroneous, being a factor of 2/ti too small. 

If the source generates isotropic radiation in the comoving frame of the expansion with a power- 
law emission spectrum hc{£c, /Uc) oc e~° , where £c and /ic are the photon energy and emission 
angle in the comoving frame, then it follows that the photon emission rate in the observer's frame 
takes the form 

n(e, /i, r) = .^£-"(1 -/3/i)-(^+") , e_ < e < e+ , (15) 

where £± define the bounds to the observed source spectrum. This form is derived from the Lorentz 
transformation relationships = re(l — [3fi) and = (// — /3)/(l — fSfi) and their associated 
Jacobian decdiid de = [r(l — /3//)]~^ , given that the total photon number n(e, /x, r)ded^dVdt 
is a Lorentz invariant. Note that in Eq. ([l5|), a factor of l/F has been absorbed in the definition 
of N . 



The value of the coefficient N in Eq. (|T5|) can be determined by computing the photon flux 
JF at large distances from the source and equating the result to the observed flux in individual 
GRBs. Specifically, the flux at test photon energy = 1 (i.e. 511 keV) is 

= cj dfiti lJ-tin{et = l, flu, n) = Af J _ ^^.-^a+l dfJ-u dm , (16) 

after integrating Eq. (|l^) over azimuthal angles. The units of J- are photons per square centimetre 
per second. For the moment assume that there are no angular restrictions to the phase space, i.e. 
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©B = T^/'2. Further, note that in this integral the test photon acts purely as a position marker 
and can be taken to be on-axis without loss of generality, i.e. 0^ = . The angle cosine /Uj in the 



distribution is given by Eq. (22) below, and the limits on the integrals are defined by the radial and 



causality constraints in Eqs. (|9|) and (10). At large distances from the source, these restrictions 
imply that ru/rt ~ 1 and //ti ~ 1 , as is obvious from the description of the geometry in Fig. |^a. 
In fact, defining ^ = rt/ Rq , then Eqs. (^) and (|lO|) can be expressed as ip- ^ "0 ^ '4'+ i where 



are given in Eq. (29). The evaluation of the integrals in Eq. (^) can be facilitated by changing 
variables thus: the integration is performed using the variable s = \/l — 2p^ti + /(I — p) for 
p = r^j/rt ( ~ 1 ), and the p integration is calculated using the variable t = 1^(1 — p) . For ip ^ 1 , 
then /ij ~ 1/s . Reversing the order of integration yields l<t<(l — /? + AR/Rq)/{s — (5) and a 
trivial result for the t integral, so that 

Here rt is set equal to the distance d between the source and the observer, and the flux naturally 
obeys an inverse-square law: oc d~'^ . Note that the algebraic manipulations here closely resemble 
those applied to the expression for the optical depth; for this reason, detail is minimized here and 
is deferred to section 2.2.2 below. 



In general, the result in Eq. ( p^ can be expressed in terms of a hypergeometric function of 
two variables, however it is simple to evaluate it directly by numerical integration. In the special 
cases where the expansion is, in the comoving frame, a thin spherical shell with AR/Rq ^ 1 — /? , 
or it is a filled shell with AR/Rq ^ 1 — /3 , it becomes analytically tractable, giving 



^ M {ARf 



— M « 1 _ ^ 



(l-/3)"+2 ' Ro 
3(a + 3)/3 l(l-/?r+^ V' Ro^ ^ 



(18) 



The quadratic dependence of J- in AR when AR/Rq <^ 1 reflects the two dimensions of the 
integration in Eq. (H). At the same time, J- is independent of Rq when Rq ^ AR since the 
inner radius contributes negligibly to the source volume. The value of M is therefore determined 



by equating the flux in Eq. (17) to that observed at 511 keV for sources with unbroken power-law 
spectra, or by a power-law extrapolation of the high energy spectrum down to 511 keV for those 
sources with spectral breaks above this energy (e.g. GRB940217; see the discussion in Section 3). 

The modification to the expression in Eq. ( [l7| ) for the flux that is induced by reduction of 
0B below tt/2 is straightforward. The considerations of angular constraints in subsection 2.1 lead 
to the simple expression for the restriction of the {pu, p) phase space in Eq. (p^). Since the 
flux is observed at infinity, pu ~ 1 and we take 0t = for the flux calculation. Then Eq. (12) 



immediately implies that s < 1/C , for ( defined just below in Eq. ([l9D, and it quickly follows that 
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sub-spherical expansions generate a flux given by Eq. (^) but with 

s+ = minjl-l--^^, i| , C = cos(^min{7r, 6b}) (19) 

substituted as the upper limit to the integral. Clearly then, the angular restrictions play no role in 
determining the flux until the solid angle [27r(l — cos 0b)] of the expansion becomes comparable 
to the fractional shell thickness AR/Rq . 



2.2.2. Optical depth for radial propagation of test photons 

The differential optical depth in Eq. (^) can now be expressed in more explicit form using 
Eq. ( p!^ and the form of the photon emission rate in Eq. (^), evaluated at the position of the 
interacting photon: 

= / (1 ! ^^^'f +1 de^ duu dru d^u ■ (20) 

Here M has been removed from the integration because it is assumed to depend only on 0b and 
be independent of the position within the source. Hereafter it will be assumed that the emission 
spectrum in Eq. (|l5|) has a large or infinite range ( e+ ^ e_ ), for which it is possible evaluate the Si 
integration separately, and analytically. Specifically, it is permissible to change variables in Eq. ( pO|) 
to the CM frame energy variable x = V^t^ii^ — IJ-u)/'^ via 4xc?x = ^t(l ~ t^ti)d£i , following the 
procedure of Gould and Schreder (1967; see also Baring, 1993), and perform the integration of the 
cross-section separately. Consequently, the differential optical depth assumes the form 

where the integration of the cross-section over x is 

nia) = ^J^ ^^-^-a,,{x)dx ^ . (22) 

The approximation in Eq. (^) was obtained (see Baring, 1993) from Eq. B6 of Svensson (1987), 
who also gave the exact analytic expression for the integral; it is accurate to better than 1% for 



1 < a < 7. The angle cosine = cos 9i of the interacting photon that appears in Eq. (21) can be 
determined explicitly from the geometry in Fig. ^3 using Eq. (0); eliminating the variables and 
Qti gives (for 9t = 0) 

^. ^ ^'^i-p . (23) 



The scaling variable p = ru jrt ( < 1 for = 1 ) will be of use in the development of the optical 
depth integration. 
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In applications where the emission spectrum is of finite energy range, the modification for 
performing the Ei integration has been developed by Gould and Schreder (1967) and Krolik and 
Pier (1991). Low energy spectral turnovers or cutoffs are unlikely to be influential in pair production 
continuum attenuation calculations applied to gamma-ray bursts that have bulk motions with large 
Lorentz factors T and maximum energies under 100 MeV (Baring, 1994). While turnovers are 
observed as photon energies drop into the BATSE range, sharp cutoffs can presumably be only in 
the soft X-ray range where photons interact with gamma-rays of energy much more than 1 GeV to 
produce pairs. Therefore any suppression of 77 e^e~ continuum attenuation that is introduced 
by a low-energy cutoff is unlikely to be observed by EGRET. However, for observations in the super- 
GeV range, spectral structure in the BATSE range becomes quite relevant to opacity determinations 
(Baring and Harding 1997), and is discussed briefly at the end of Section 3. Introduction of high 
energy cutoffs to the spectrum of interacting photons is also largely irrelevant to this investigation 
because the most energetic EGRET source photons predominantly interact with photons at energies 
considerably below the maximum detected. 

The azimuthal integration in Eq. (^) can be performed after establishing the restrictions 
the source geometry places on (pu ■ For the interacting photon to be in the expanding "conical 
shell," the radial restriction in Eq. (P) is independent of (f>ti ■ In contrast, the angular constraint 



< Qi < @B in Eq. (11), which is independent of time, does restrict the allowable azimuthal 
angles. Assuming that the plasma emits uniformly at any one time, this angular constraint results 
in an analytic determination of the azimuthal integration, because at a given radius each azimuthal 
angle within the cone of expansion contributes equally. Since the azimuthal angle (pu is defined 
to be the angle between the planes OZT and OTI (see Fig. |l|b), then considerations of spherical 
trigonometry yield 

cos Qi = cos Qt cos Qti + sin 6^ sin Qu cos (pa > cos 0b (24) 

when |0tib0tj| < 0B . When this condition is not satisfied, the interacting photon is always within 
the cone defined by the expansion and all values < < 2tt are permitted. It follows that the 
(pti integration has limits defined by | cos (pul < , where 



COS 9b - cos Qt cos Qti p> r, ^ n ^ n ^ r\ 

sin sin G*. ' Ob - < O,. < Ob + , ^^^^ 



1 , 0<Gti<GB-Gt 



The restriction Qa > Qb — Qt is necessary to achieve 77 > —1, and as Qa — > Gb -|- Gt then 
7/ — > 1 and the permitted cpti phase space shrinks to zero. When Gf — > , (i.e. the test photon 
is on the axis of the expansion), r] ^ —1 for all permissible Qa ■ This simple special case will 
be assumed throughout subsequent sections of the paper. The (pti integration in Eq. (^) is then 
trivially evaluated to give 2 arccos r] . Equation (^Tj) can then be integrated over the test photon 
position rt to give the total optical depth: 

T,,{et) = ^ ^ er' ^ J arccos r, ^[^J^^^f+i df,t, dn, dn , (26) 
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and the value for t] when \Qt i @ti\ < ©b becomes (for 9t = ] 



cos9bV1 - 2p/itj + - cos Gt(l - p/Uii) 



psinGtA/l-/ii^i 



where the substitution for @ti has been effected using Eq. (^), the variable p is defined in Eq. (13), 
and the sine rule applied to triangle OTI in Fig. |l]b. 

For the moment, consider 47r steradian expansions, where 0b = vr and angular restrictions to 
the interaction phase space do not enter the analysis. Evaluation of the triple integral in Eq. ( p6D 
can be facilitated by changing variables via the scaling transformations defined in Eq. ([l^ ) , so that 
Hi and T] are rendered independent of rt (see Eqs. (p3|) and (^)). A complete set of dimensionless 
variables has therefore been chosen and reversing the order of integration so that the ip integration 
is performed first yields an analytic reduction of the optical depth to a double integral: 

Tyyiet) = M^Rhr' ^ fdp f dpu {^l-i^D arccosry-i^^^ . (28) 



The value of //min is given in Eq. ( pOj ) . Here il)± are the limits of the ^ integration ( < ip < ip+ ) , 
and are determined from the radial constraint in Eq. (^) and the causality condition in Eq. (p!0|): 



where < /) < 1 defines a "look-back" volume. If the expansion is at least fully hemispherical 
( O > 27r ) , then rj = —1 , and the pti integration is over a range determined by the condition 



-(/;+ > ip- ■ It follows from Eq. (|29| ) that the range of pu permitted in Eq. (28) is Pmin ^ Pu ^ ^ , 
where, for any P < 1 , 

.„„ ^^{'^''■^-"-'■^^^fn- M<^^'- (30) 

1—1 , otherwise. 

Clearly Pmin < 1 is always true. This restriction indicates that the look-back volume is not 
spherical because of the presence of edges to the expansion in the radial direction. The form of the 
pair production optical depth in Eq. ( p8D is not generally reducible to a simpler form, and is ready 
for numerical evaluation in cases where 0b > vr/2. However, analytic development is possible in 
the special case of propagation of test photons along the axis of the expansion (i.e. 0t = ), which 
will be treated in the next subsection. 

Consider now the additional restrictions on the integration phase space due to a reduction in 
the expansion opening angle ©b • The way the expansion has been defined automatically precludes 
any necessity to treat cases where 0b > 7r/2 , since they reduce to the 6b = n/2 case. This 
arises because switching on the expansion at t = implies that interacting photons from the back 
hemisphere 0j > 7r/2 (for 0t = ) can never reach the test photon originating in the forward 
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hemisphere, and therefore cannot contribute to the optical depth. This switch-on stipulation is a 
reasonable approximation to a real burst, and the contribution to the optical depth from photons 
originating in the back hemisphere is expected to be strongly suppressed due to the relativistic 
nature of the expansion. When 0b < vr , both nu and the azimuthal angle cpu are restricted. The 



constraint on Qu in Eq. (12) can be inverted to find the ranges of acceptable values of iJ-u in terms 
of p, which defines how the reduction of the opening angle of the expansion restricts the {f^u, p) 
phase space: 

- 1 < Aiti < ; /"+ < ^ti < 1 , (31a) 



where 



M± = <! ^ ^ ^ (31b) 



Therefore when p < ^1 — all values of pti are permitted, since then the angular boundary 
to the expansion lies outside the look-back volume defined by p . The azimuthal restrictions are 
reflected in the value of rj in Eq. (p7[), and by a similar analysis, the boundary where r] increases 



above —1 is also defined by Eq. (|3lD , but with ( cos(Gb — @t) ■ The ranges in Eq. (|31|) are, of 



course, subject to the Pmin < Pu limitation that is imposed by radial considerations. 

In the case where the test photon propagates radially along the axis of the expansion (i.e. 
Qt = 0), which will be the focal point of all subsequent developments in this paper, the double 
integral expression for the optical depth in Eq. (^) can be manipulated into a form that is more 
convenient for numerical computation, where one integral can be expressed in terms of familiar 
hypergeometric functions. The analytical approach is similar to the derivation of the photon flux 
in Eq. (|l^) . The pti integration is expressed in terms of the variable s , defined in Eq. ( |l^ ) , and 
the p integration can be performed first using the variable t = 1 — p . This change of variables 
leads to the range 

A 7? 

1 < s < 1 + ^ , (32) 
no 

for the variable s , which can be easily deduced from the requirement that pti < 1 and the condition 
that '(/'+ > "0- ill Eq. (|29|). For specific s within this range, inversion of the restriction in Eq. ( ^0|) 
leads to the upper bound t+ = 2/(1 -|- s) for t, which can also be obtained equivalently from 
the radial constraint in Eq. (^). The lower bound to t can be derived directly from the angular 
constraint in Eq. (12); the result gives 



max|o,^^-^| = t_ < t < i+ = ^ , (33) 

where C is defined in Eq. (^). Since C ^ 1 1 it follows that t+ > t_ so that only one range arises 
for the t -integration. This simplicity does not arise if the s -integration is performed first since 
t_ has a maximum of 1 — v^l — (^^ that can then yield two integration ranges for s for some t . 
Note that t_ exceeds zero only when s > 1/C • Hence, the angular constraint only impacts the 
calculation of the optical depth when > Rq/{Rq + Ai?) , a situation identical to that arising in 
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the treatment of the source flux. Remembering that arccos rj = tt for this case of radial and axial 
propagation of test photons, the optical depth in Eq. (28) develops to the form 



fjf CTt r>2 ^(«) 



l+<5 



dt 



t 



ds 

2a 



V 



T+a 



1-13 + 5 



, V = 2{s- f3)+t[il + s^)p-2s] 



(34) 



for fractional shell thickness 5 = AR/Rq 



The angular constraint in Eq. (|12|) is unimportant only when t_ = . In general, this is not 
so, and the t -integration in Eq. ( p^ can be written as the difference between integrations over the 
ranges [0, t+] and [0, t_] . Therefore two terms appear, each of which can be manipulated in a 
similar fashion. Consider first the integration over [0, t_] : the t -integration can be rewritten via 
the substitution t = t-{l — q) , leading to the transformation 



{s^ - ly 



dt 



t 



2a 



2"cj"+i S 
\«+r 



Here A = A(s, /3) and o = (7(s, /?, C) are given by 



a A 



dq 



(1 - aAg)"+^ 



, _ (l + g^)/3-2g 
^ - s{s-l){l + (3) 



a 



(s-l)-(l-C)sA 



(35) 



(36) 



It can be shown that A is a monotonically decreasing function of s with the range — oo < A < 
/3/(l -I- P) , and furthermore that < a < 1 . When C ^ 1 (i-e- ©b ^0), then and 
cr ^ 1 , so that the result for t-integration over the range [0, t+] is recovered. Putting the two 
terms together and defining the Heaviside step function TCs{x) to be unity when x > and zero 
otherwise, the optical depth can therefore be written in the form 



_ _ jCr CTt ry2 a-l _H(a} 



l+<5 



ds ■ 



I- [3 + 5 
s- 13 



J„(s; /3, C) (37a) 



for 5 = AR/Rq , with 

Ms; PA) ^ aa(A)-Ws(s-i) (^^Y a'^^^ g^{cj\) 



(37b) 



The function Qa{z) is just the integral that appears in Eq. (35), and is expressible in terms of 
the standard hyper geometric function F{a, (3; 7; z) using the identity 3.197.3 in Gradshteyn and 
Ryzhik (1980): 



Qa{z) 



dq 



[I 



\2a 



1 



(1 - zq) 



a+T 



l + 2a 



F{a + 1, 1; 2a + 2; z) 



(38) 



The numerical evaluation of Qa is straightforward, and is described in the Appendix. An alternative 
form for the optical depth can be derived by leaving the t-integration in Eq. (34) as one integral 
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over the range t+] and rescaling the integration variable. This second form yields Eq. (|37|a) 
but with an alternative representation of the function Ja{s] f3, C) '■ 

Ms; P, C) = /^^ dq , g+ = minjl, ^^^1 . (39) 

While slightly less convenient than Eq. ( p7| ) for numerical evaluation of the optical depth, this 
second form is useful when obtaining results in the limiting case of small opening angles, treated 
in subsection 2.3.3 below. Computationally, if the series representation of Ga described in the 
Appendix is used, Eq. (|37|) is a single integral that is relatively simple to evaluate. Note that the 
integrand does not diverge at s = 1 due to the behaviour of GaW there (see the Appendix). 

It is important to emphasize that the optical depth in Eq. (^) was obtained under the assump- 
tions that the test photon originates at the rear of the expansion and propagates radially outward. 
With this specification, Eq. (|37|) is intended to approximate a variety of possibilities for test photon 
initial conditions. Non-radial test photon motion will increase the optical above that in Eq. (^7|), 
primarily because of increased angles with interacting photons, however an observer's unique per- 
spective will strongly bias against such situations for relativistic expansions. On the other hand, 
test photons can plausibly originate closer to the surface of the expansion than Rq , diminishing 
the optical depth accordingly. Suppose that the test photon starts at radius Rq + I'AR at time 
t = . Then all interacting photons inside this radius are always causally-disconnected from the 
test photon because it propagates radially. Hence the region interior to Rq + i^AR is irrelevant 
to the optical depth calculation, and a new initial condition can be defined, with Rq + I'AR and 
Rq + AR denoting the relevant inner and outer radii of the expansion. The test photon is now 
at the rear of this section of the conical shell, and the optical depth computation can be repeated 
entirely with the aid of the substitution 

Rq Rq + vAR, AR {l-u)AR (40) 

without any additional manipulation. This elementary transformation propagates all the way 
through the development so that the optical depth for test photons starting at arbitrary positions 
Rq + uAR in the expansion is simply from Eq. (^) by the substitution 5 — > (1 — v)5/{l + v5) , a 
very attractive scheme of generalization. Note that the normalizing flux in Eq. (|l^) is unaffected 
by these considerations. 

This concludes the analytic development of the optical depth formula, which is ready for 
numerical computation and certainly is much more amenable than the quintuple integral in Eq. (^0|) . 
Before presenting such computations (in subsection 2.4 below), it is instructive to examine the 
optical depth for some limiting cases of the source geometry. 
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2.3. Approximations in Limiting Cases 

There are four special cases where it is both possible and useful to obtain analytic limits to the 
pair production optical depth: these correspond to the thin-shell limit, thick-shell or filled-sphere 
expansions, narrow beams, and a stationary photon gas. 



2.3.1. The thin-shell limit 

The expression in Eq. (|37| ) is in suitable form for the derivation of the optical depth in certain 
special cases. The first of these is the limit of a thin, spherical shell for the expansion, where 
5 = AR/Rq <C 1 — /3 (i.e. the shell is also thin in the comoving frame of the expansion), and 



(5 <C 1 — C so that the angular constraints are immaterial. The ^a(<TA) term in Eq. (37) is 
therefore absent. As noted in the Appendix, in the s — > 1 limit, (1 + 2a)^o(A) approaches 
F{a + 1, 1; 2a + 2; !)/(! — A) , which with the aid of identity 9.122 of Gradshteyn and Ryzhik 
(1980) leads to the limit ^^(A) l/[a(l - A)] « (s - 1)(1 + (3)/[2a{l - /?)] . It is then elementary 
to derive the result T^^{et) oc M {AR^ . The two powers of AR are due to the thin shell severely 
restricting the spatial extent of the and integrations [e.g. see Eq. (]20|)]. Such behaviour is 
largely meaningless until the dependence of the formula for the flux is factored in. Remembering 
that in this limit, for a fixed observed fiux, Eq. ([l8|) yields M oc l/(Ai?)^ , the optical depth is 
virtually independent of AR , as is expected. Explicitly, we obtain 

r..(e.) - , ^«min{l-/?,l-C} . (41) 



This thin-shell limit displays a strong inverse dependence on the bulk Lorentz factor F = — 
of the expansion, as is evident from Fig. ^: typically a ~ 2 — 3 for EGRET bursts (e.g. see Table |2|). 
This dependence is enhanced by one or two powers of F that appear in the determination of Rq . 
This case is closest to the work of Fenimore, Epstein and Ho (1993), who treat test photons coming 
from an entire spherical shell (i.e. including the limb regions). As argued in Section 2.1, the major 
contribution to the optical depth comes from test photons originating in near-axis environs so that 
the differences between conclusions made using a formula like Eq. (^Tj) and the work of Fenimore 
et al. (1993) are marginal. 



2.3.2. Filled-sphere expansions 

The other extreme class of expansions from the point of view of the radial dimension contains 
filled spheres initially, i.e. thick shells (in the comoving frame) with AR/Rq ^1-/3. Again we 
shall ignore the impact of narrowing the solid angle down and demand 0b = 7r/2 here, so that 
only one term in Eq. ( ^7|) contributes to the optical depth. In this limit, inspection of Eq. ( ^7|) 
soon reveals that the dominant contribution to the integral is for s — 1 <^ \ — (3 . Then it follows 



- 19 - 



that A ~ — (1 — — 1), since relativistic expansions with /? ~ 1 are considered here. Choosing 
A/(A — 1) as the integration variable and using the transformation of the hypergeometric function 
in Eq. ( |A2| ) yields a result that is proportional to the integral in Eq. ( |A4| ). As in the thin-shell 
limit, here r-y-y(ej) oc (Ai?)^ , reflecting the dimensionality of the integrations. The filled sphere 
limit of Eq. (^) indicates that an observed flux gives a volume-determined photon injection rate 
M oc l/(Ai?)3 , so that the overall expression for the optical depth is a declining function of the 
expansion thickness: for /? ~ 1 , the asymptotic result 



'77 V 



i+1 ^a-l 



(a + 3) 



a a 



{^(2a) - ^(a) - l} 



(42) 



is derived, where 'ip{x) is the derivative of the logarithm of the Gamma function, defined in Eq. ( A4). 
Again, a strong dependence on the bulk Lorentz factor of the expansion is evident. The moderate 
decline of T^-y{£t) with Aii reflects the fact that large regions are less compact for a given source 
luminosity. Note that extremely-filled spheres with Ai? > Rq are not really discussed in this paper; 
these seem unlikely to be realized in bursts and require an alternative coupling of length scale to 
time variability, via Ai2 = cAi , as is mentioned below. 



2.3.3. Narrow beam expansions 

The case of small opening angles of the expansion is also of interest. Since the emphasis here 
is on axial viewing perspective, this limiting case corresponds to ^ = cos Gb ~ 1 . Specifically, 
this narrow beam limit satisfies 1 — C<Cl — /3,so that the reduction in opening angle dominates 
the causality limitations, and also 1 — C, 6 = AR/Rq. However, such an identification with 
small solid-angles is not sufficient to define narrow beam cases; as will be evident shortly, the 
size of the opening angle itself is also quite pertinent. The most suitable form of the optical 
depth for development here is using Eq. (37 a) combined with the representation in Eq. (p9|). The 



s -integration then consists of a range 1 < s < 1/C over which the volume is not opening angle- 
limited, and this contributes of the order of 1 — C to the optical depth expression in Eq. (|37|a). 
However, it turns out that the range 1 — C s — 1/C ^ 1 dominates the contribution to the integral. 
This leads to a simplification for the integral in Eq. (|39|) for J^a{s; /?, C) • Yet the integral is not 
trivial since the parameter Ag+ is not necessarily small; in fact —Xq+ ~ s(l — C)(l — — 1)^ ■ 

Reversing the order of integration and using 3.194 of Gradshteyn and Ryzhik (1980) leads to the 
result that T^^{et) oc M^/l — ( . This result is subject to the requirements that vT^^ <^ 1 — f] 
and ^yl — ( <^ 6 . Clearly the reduction in angular phase space for pair creation is not solid-angle 
limited, but rather constrained purely by the maximum size of the angle between the momentum 
vectors of the test and interacting photons — this scales as 0b . Such a linear dependence on 
Gb is manifested in both this asymptotic limit for Eq. (|37|a) and the range of its validity. When 
6 <C VI ~ C then the optical depth in Eq. (^a) no longer behaves like a narrow beam limit, but 
rather like a thin-shell limit as discussed above, where r^-y 
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The other part of the narrow beam calculation relates to the flux normalization factor T 
in Eq. ([l^). Remembering the limitation in Eq. (p!9D, the computation of J- is straightforward, 
yielding a solid-angle limited flux: T oc AA(1 — Q . This intuitively obvious result is applicable only 
when 1 — C ^ 1 ~ /3 and 1 — C <C 5 . Such a solid-angle-limited range differs from the requirements 
imposed by Eq. (37), thereby complicating the consideration of the narrow beam limit. To aid 



understanding of this limit, the various dependences of Eq. (|17|) and Eq. (p7|), and the resulting 
behaviour of the overall optical depth, as functions of 0b and b , are listed in Table |l|. There, four 
parameter regimes with \/\ — Q ^1-/3 are identified, depending on b . Three of these regimes 
are strictly narrow beam limits, while the fourth, for 5 <C 1 — C j corresponds to the thin-shell limit 
in Eq. (41), and is independent of 0b . For the first two regimes in Table ||, the developments of 



Eqs. ( ITTD and ( |37| ) just mentioned lead to an approximate overall optical depth that can be written 
as one expression: 

72 



, . -ia^d^^ r- r(a + 1/2) n{a) [2(1 - + 5] 

- 2c Ro-^"^^ r(a + l) 3(1-/?)^ + 3(1-/3)5 + 5^ 



for 5 = AR/Rq . This formula encompasses both thin-shell, narrow beam ( ^/T^^ <C <C 1 — /3 ) 
and thick-shell, narrow beam ( y/1 — C <C 1 — /? <C (5 ) regimes. Surprisingly, the calculated optical 
depth actually increases when the opening angle closes down, reflecting the explicit dependence [see 
Eqs. ( p!?! ) and ([l9|) ] of the observed flux on the solid-angle of the expansion, combined with the pair 
production rate explicitly depending only on the angle between the test and interacting photon 
momenta. Essentially, the photon density in the source increases for constant observed flux as the 
expansion opening angle is reduced. The resulting optical depth varies only with the thickness of 
the shell (oc 1/AR) in the regime of thick shell expansions, i.e. for AR/Rq 1 — /3 . The third 
regime in Table [l| requires use of the thin-shell evaluation of T^^jM , and yields the asymptotic 
approximation 

r,,(.0-2x|.^,«-.«M('l_|)"_A_ , l_c«i«v^«l-,3. (44) 



This bears an even stronger dependence on the opening angle Gb ~ \/2(l — ^) , again rising with 
increased narrowness of the beam, a property that corresponds to an enhanced mean density of 
radiation in the expansion. In this case, the phase-space for pair production is not restricted by 
the opening angle, but rather only by the thinness of the shell, while the flux is still solid-angle 
limited. 

Finally, in concluding the consideration of narrow beam cases, observe that Eq. (^) approx- 
imately reproduces the thin shell and filled-sphere forms in Eqs. (|4l| ) and ( |4^ ) when the beam is 
opened up to 1 — ^~1 — /?. For this intermediate (or critical) regime of opening angles, a domain 
common to all three of the limiting cases discussed so far is achieved when AR/Rq ~ 1 — /? . This 
locality in phase space corresponds to the so-called "blob" scenario of earlier work (e.g. Krolik and 
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Pier 1991; Baring 1993; Baring and Harding 1993) on pair production transparency constraints in 
gamma-ray bursts, a situation that is discussed in Section 3. 



2.3.4- Stationary radiation gas 



The remaining limiting case of the optical depth is for a /? = or non-relativistic expansion. 
This is mostly of academic interest as a check on the numerical evaluations that follow, and does not 
have great physical import for the problem considered in this paper. In fact the result in Eq. ( ^7\ j 
does not aptly model stationary gases since we have neglected the limb contributions to the optical 
depth; these become significant in non-relativistic expansions. Set ©b = vr/2 for simplicity. In the 
(5^0 limit, A — 2/(s — 1) , which simplifies Eq. (]37| ) somewhat. However, analytic development 
of the subsequent result is not fruitful, so it is convenient to consider separately the thin and thick 
shell cases. The limiting result in Eq. (^) was derived without restriction on (3 , so the limit /? — > 
can be taken to obtain the optical depth for stationary, thin shell sources. For AR/Rq 3> 1 — 
a derivation alternative to that in subsection 2.3.2 is requisite. Then the dominant contribution 
to the integration comes from s — 1 <^ 6 . The transformation relation in Eq. ( [A2| ) can be used, 
together with a change of variable to yield a result proportional to both AR'^ and the integral in 
Eq. (|A5|) . The flux is simply evaluated when /3 = , so that the optical depth for a thick stationary 
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Note that this estimate differs somewhat from results derived for isotropic photons (e.g. Schmidt 
1978; Epstein 1985; in particular those that use the formalism of Gould and Schreder 1967), because 
in this paper we have assumed isotropic injection in the comoving frame (in this particular limit 
the observer's frame), which is not equivalent to radiation isotropy due to radiative transfer in the 
sphere. 



2.4. Numerical computation of the optical depth 

The various limiting cases just explored guide the technique for numerical determination of the 
optical depth, and further act as checks on computational accuracy. The numerics are generally 
straightforward, and it is expedient to use s — 1 as an integration variable in Eqs. (|l7|) and (^7|), 
and scale quantities in terms of 1 — /? to maintain good accuracy for large Lorentz factors. Since the 
range of integration variable contributing significantly to the two integrals is sometimes quite large, 
logarithmic sampling of s — 1 is favoured. Use of the functional form in Eq. (|37|b) comfortably 
produces smooth results for ©b down to even smaller than 0.01 degrees, so use of the alternative 
representation in Eq. (39), or a series in 1 — , is unnecessary for the purposes of this paper. 
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Numerical determinations of the optical depth formed by the combination of Eqs. (^) and 
( [l7| ) are presented in Fig. |2| for different fractional shell thicknesses AR/Rq illustrating its strong 
dependence on T , and significant dependence on source opening angle only for smaller T . The 
inclusion of /? ^ 1 cases is intended only to provide a general guide to the behaviour of the optical 
depth (perhaps to order of magnitude accuracy), but this is strictly incorrect since they neglect 
limb contributions. The quantity actually plotted in Fig. |2| is the scaled optical depth 



T-y-y(et) At (ms) 



= 1.06 X 10^3 ^ 7^(a; /3, AR/Rq) , (46) 

/in 



which is dimensionless since J- is measured in photons per square centimeter per second, where 
7^(a; /3, AR/Rq) is just 7i(a;)/(l + /5)i+° times the ratio of the two integrals in Eqs. (|7|) and (|l^. 
Here At is the observed source variability timescale, typically in the range of 10~^ -1 seconds, 
inferred, for example, from time histories such as those exhibited in the BATSE IB catalogue 
(Fishman et al. 1994), and z is the cosmological redshift of the source (it is set to zero in Fig. |2|). 
The coefficient of this equation clearly defines the optical depth scale for cosmological bursts, and 
would be 8-10 orders of magnitude smaller for galactic halo sources. A "canonical" spectral index 
of a = 2.5 (see Table |2| for specific values) is chosen in Fig. ^ for simplicity; increasing (reducing) 
a just increases (lowers) the slope of the curves in the F S> 1 regimes. Remembering that the 
energy £t is expressed in units of mgC^ , it is evident from Fig. ^ that the maximum observed 
energies of EGRET sources (see Table ^) imply that Lorentz factors in the range of approximately 
50-500 are required to render these bursts optical thin to pair production. 

The curves in Fig. ^a clearly delineate three regimes of parameter space in order of increasing 
bulk Lorentz factor: (i) non-relativistic flows where the optical depth is independent of the expan- 
sion speed, (ii) thin-shell expansions ( AR/Rq <^ 1 — (3 and F ^ 1 ), which is the portion of phase 
space sampled by the work of Fenimore, Epstein and Ho (1993), yielding a strong inverse dependence 
of r^-y on F , and (iii) at the highest Lorentz factors (i.e. above the break at (Ai?/i?o)/(l — /3) ~ 1 ), 
thick-shell expansions, where the inverse dependence on F is slightly stronger (by two powers of 
F : compare Eqs. (||) where r^^ oc (1 - /?)"/i?o and (|4|) where r^^ oc (1 - /?)"+VAi2). We note 
that Fenimore, Epstein & Ho (1993) produced an optical depth vs. logF plot for their infinitely 
thin shell analysis that exhibited a dramatic reduction of r-y-y above some critical Lorentz factor. 
These turnovers were found to be artificial, being caused by a coding error (Fenimore, private 
communication). The ratio AR/Rq is independent of F in Fig. though other choices are quite 
plausible. For opening angles 0b ^ l/F, the ratio x = (AR/Rq) /{I — /?) is the only critical pa- 
rameter delineating the thin-shell and thick-shell cases, and is intimately related to the the portion 
of the emission region that is causally-connected to the test photon. As is evident from Eqs. ( p|) 
and (^), the optical depth is independent of the thickness of the shell when this parameter is 
much less than unity, and inversely proportional to AR when x ^ 1 i the transition between 
these regimes is quite gradual. The curves exhibit a lack of dependence on the opening angle 0b 
when 0B ^ min{l/F, {AR/ Rq)^/"^ } (the curves in Fig. ^ are coincident for 0b = 90°, 10° , and 
1°): this is a consequence of causality dominating geometry in the restriction of the integration 
phase space. When 0b ^ 1/r, the optical depth (oc (1 — /3)""'"^/^/[0b -Ro] ) is actually increased 
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by the angular reduction of the available phase space, an effect that is particularly evident for 
non-relativistic expansions. As mentioned in Section 2.3.3, the reason for this increase in optical 
depth with declining is that the pair production rate is proportional to the angle between the 
test and interacting photons (which is limited linearly by Gb), while the flux scales as the solid 
angle (i.e. 0^ when @b <C 1 ) of the expansion. Hence the density of photons in the source that is 
inferred for a given flux actually increases as 1/0b when Bb declines below 1/r. In Fig. ^3, where 
AR/Rq is approaching unity, the thin-shell regime is no longer distinct from the non-relativistic 
domain so that all relativistic expansions are thick-shell. 

The specific choice of Rq = TcAt is made in Fig. § for simplicity; alternative dependences 
on r are possible, corresponding to different interpretations of source timescales (as discussed in 
Section 3 below), and these result in only a slightly different appearance from the curves in Fig. |^. 
The results presented in this subsection are obtained under the assumption that the test photon 
starts its life at the rear of the expansion; Eq. ( |40| ) provides a simple scheme of substitution for 
Eq. ( [37| ) that yields optical depths for arbitrary initial positions of the test photon within the 
expansion. It is evident that imposing a T-y-y = 1 condition on EGRET bursts potentially can map 
over into both thin and thick shell regimes, depending on the assumed source distance, measured 
flux and maximum photon energy observed. We note that using a r-y^ = 1 criterion for source 
transparency actually leads to conservative lower bounds for F , since significant attenuation is 
already present when r^^ = 1 . This can be easily seen if the spectrum is attenuated by an 
exponential factor exp{— r^^} , a common choice. A central consequence of the assumption of 
an infinite power-law burst spectrum is that the optical depth is an increasing (and power-law) 
function of energy, so that spectral attenuation arises only at high energies; departures from this 
behaviour will be discussed briefiy at the end of the next Section. We remark also that opacity skin 
effects (which depend on the spatial distribution of photons and therefore are model-dependent) can 
sometimes render the exponential exp{— r^^} a poor descriptor of attenuation, with 1/(1 -|- r^^) 
perhaps being an improvement for uniformly-distributed photons, leading to broken power-laws 
rather than exponential turnovers. A variety of signatures of spectral attenuation are possible, 
particularly if pair cascading is involved, and some of these are illustrated in the work of Baring 
and Harding (1997). 

Undoubtedly the most crucial piece of information to be gleaned from Fig. |2| is that causality 
minimizes the role of the opening angle 0b of the expansion in faster expansions, and specifically 
that the r^^ = 1 condition will be virtually independent of 0b in the range 90° > 0b ^ l/F 
for burst sources. This insensitivity to the angular extent of the expansion is the keystone to 
elimination of the number problem for cosmological bursts, as is discussed below. 
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3. BULK LORENTZ FACTORS FOR EGRET SOURCES 

The calculations for r^^ we have performed can readily by applied to gamma-ray bursts 



detected by EGRET. The pair production optical depth in Eq. (37) depends on the free parameters: 



r, AR/Rq, 0b, d, and on the observed parameters: source flux [defined in Eq. (|T^)], high 



energy spectral index a [see Eq. (15)], test photon energy £t (= Emax), and At. The burst 
variability timescale At can be used to infer an upper limit on the source size, which we take for 
the moment to be Rq = TcAt , based on the apparent size of the expanding shell perpendicular to 
the light of sight, as seen by a stationary observer (e.g. see Rees, 1966). Alternative size-estimates, 
such as those that couple the source variability to dimensions along the line of sight to the source, are 
possible. In addition, measured minimum variability timescales for bursts range from milliseconds 
in the BATSE range to supersecond values in EGRET data. Motivations for choosing either of these 
At , and also different source size determinations, are discussed below. Both variability timescales 
are addressed in the results presented here, for the sake of completeness. 

The key observable parameters, to be used in the pair production opacity calculations of this 
paper, are displayed in Table § for six of the burst detections by EGRET. Note that there are 
eleven EGRET bursts in total (Schneid et al. 1996), three of which have insufficient (published) 
data for the purposes of our analyses; GRB 920622 and GRB 940301 have the required observational 
parameters in Schneid et al. (1995), but suffer from poor statistics above about 2 MeV. We therefore 
conservatively opt to study just the most significant six sources of the EGRET population. In 
Table |2|, the fluxes are expressed in observer- friendly units, via /(I MeV) , which is just the flux, 
evaluated at 1 MeV, per MeV energy interval (the flux is per mgC^ at 511 keV). These fluxes are 
obtained via extrapolations down to 1 MeV of the best-fit power-laws to the time-integrated super- 
MeV spectral data, and are not necessarily the actual fluxes measured at 1 MeV. Note that both 
EGRET and COMPTEL parameters are listed for GRB 910601 since this burst was relatively soft 
and actually had a slightly more significant detection by COMPTEL than by the EGRET TASC in 
the 3-5 MeV range. We remark that the COMPTEL listings for GRB 910601 (based on Hanlon et 
al. 1994) differ slightly from those quoted by Baring and Harding (1993) that were obtained from 
Winkler et al. (1993). As noted in the Table caption, in the computations of this paper we neglect 
the highest energy (18 GeV) photon detected for GRB 940217. This conservative step is taken 
because the statistically-limited sample provided by a single photon leads to a large uncertainty in 
the spectral form at these energies, which is compounded by the lack of contemporaneous spectral 
information at lower (i.e. sub-MeV) energies (e.g. Hurley et al. 1994); time-resolved spectra of 
good statistical quality in the super-100 MeV range await future generations of instrumentation. 
A nice depiction of the relative fluxes and spectra of four of these bursts is given in Hurley (1996). 
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3.1. Geometries with 0b ~ l/F 

Before presenting the results for the bulk Lorentz factor constraints inferred from our optical 
depth calculations here, it is instructive to first review the results of previous, more primitive 
bulk motion determinations. As outlined in the Introduction, the earlier work of Krolik and Pier 
(1991) and Baring (1993), and subsequent papers, considered "blobs" of radiation-emitting material 
moving at relativistic speeds more-or-less along the line of sight to an observer. The angular extent 
of these blobs, i.e. the width of the angular distribution of photons as measured in the observer's 
frame, was assumed to be comparable to l/F , where T was the Lorentz factor of the blob ( F ^ 1 ). 
For such relativistically-moving blobs, the minimum bulk Lorentz factors at redshift z are obtained 
(i.e. for T^^ = 1 ) from the result derived by Baring and Harding (1993, corrected in Harding 1994): 

rl+2° > (3.83)" (1 + z)""^ ^kpc ( gMAX f(-, yv y-^ ( aj) 

^ 3a5/3(4/3 + a)2^/ii At(ms) llMevJ MeVj . (47) 

Here /(I MeV) is just the source flux evaluated at 1 MeV, per MeV energy interval, and therefore 
is proportional to J- . This formula can be approximately reproduced from Eq. ( p|) by setting 
Rq ~ TcAt . Alternatively the thick shell approximation in Eq. ( ^2|) gives more or less the same 
estimate if AR ~ cAt/F . It therefore follows that the "blob calculation" corresponds to the 
boundary between thin and thick shell cases where the line-of-sight and transverse (e.g. variability) 
timescales are comparable. Furthermore, it also coincides with the boundary of narrow beam 
expansions, namely for 0b ~ 1/F , as can be established by setting 1 — ( = 1 — P in Eq. (^) and 
choosing intermediate shell thicknesses. 

We remark here that the choice of variability timescale for use in the estimation of minimum 
source bulk Lorentz factors is subjective. The smallest At observed in BATSE data (see Fishman et 
al. 1994) is of the order of milliseconds, and it is quite conceivable that intrinsic source variability 
occurs on even shorter timescales. Such small At were adopted in the work of Baring (1993), 
Baring and Harding (1993, 1995, 1996) and Harding (1994), and lead to light-crossing time size 
determinations of cAt ~ 3 x 10^ cm. However, variability in the hard gamma-ray band, i.e. for 
Comptel and EGRET data, can only be conclusively inferred from more severely photon-limited 
samples on the order of 0.1 second to 1 second timescales. Hence a conservative approach, adopted 
for example by Ryan et al. (1994, for GRB 930131) and Winkler et al. (1995, for ORB 940217), 
uses these longer At values in obtaining pair production constraints on bulk motion in bursts. 
Note that the relative timescales in the different GRB energy bands may not be related at all to 
intrinsic source properties but merely reflect current instrumental limitations in the time domain. 
While experimentalists might prefer the conservative variability values, theorists are often motivated 
to attribute the shortest timescales to regions emitting the highest energies of radiation. This is 
frequently justifiable in astrophysics, since a whole host of cosmic objects are powered from a central 
region and thereby generate their most energetic photons closer to the center; such photons would 
generally be expected to couple to shorter timescales. If a central "powerhouse" is indeed responsible 
for gamma-ray burst emission (as suggested by Fenimore, Madras and Nayakshin 1996), then 
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sub-millisecond timescales might accurately reflect the source conditions appropriate to EGRET 
photons. In contrast, if bursts approximate more closely the class of fireball models that generate 
the emission we see by impact on the interstellar medium (e.g. Rees and Meszaros 1992, Meszaros 
and Rees 1993), perhaps through diffusive acceleration of particles at shocks, then the highest energy 
photons are produced by particles diffusing on the largest scales of the system, and therefore might 
be expected to have longer variability timescales than photons in the BATSE energy range. To 
accommodate a variety of perspectives, this paper considers At values of 1 ms and 1 second. 

The solutions of Eq. (^) for the minimum bulk Lorentz factor F^in for the burst parameters 
of Table | are listed in Table for the two different variability timescales, and for four different 
GRB source distances that are typical of galactic disk, galactic halo, and nearby and distant 
cosmological populations (though the redshift z = is chosen for simplicity). The immediately 
obvious conclusion is that, except for a disk origin of bursts, relativistic bulk motion is generally 
inferred for the EGRET sources, because of the detection of energetic photons. The Lorentz factors 
obtained for cosmological distance scales far exceed those inferred for extragalactic jets in active 
galaxies. Conversely, Table |3| indicates that isotropic emission cannot be supported in GRBs unless 
they are quite local, i.e. well within the galactic disk, a conclusion that differs from Schmidt's 
(1978) early work principally because of the positive EGRET detections in the CGRO era. These 
estimates for the minimum T in bursts are a good first guide to constraints on bulk motion in their 
emission regions; the refinements of the burst geometries addressed in this paper only modify the 
estimates in Table ^ by factors of at most a few, as will become evident shortly. We remark that 
previous versions of these estimates (e.g. Baring 1993, Baring and Harding 1993, Harding 1994, 
Baring 1995) have sometimes used slightly different observational parameters. Note that a number 
of entries in the 1 kpc column have Tmin = 1 • These actually represent unphysical T < 1 solutions 
to Eq. ( |47| ) that are obtained only because the assumption F ^> 1 is used to derive Eq. (|47|). Hence 
Fmin = 1 entries denote regimes where this "blob" constraint breaks down and refinements are 
needed. It must be emphasized that Eq. ( J^7\) and the work of this paper implicitly assume that 
the GRB spectrum extends above 511 keV in the rest frame of the emission region, so that phase 
space above pair creation threshold is non-zero. 

It is appropriate to remark on a caveat to these results. Given emission observed out to Smax , 
the maximum photon energy in the source rest frame is of the order of Smax/F , which must exceed 
unity in order to be above pair threshold. Hence the possibility that intrinsic cutoffs could be 
present in the GRB spectrum anywhere above Emax automatically implies that Emax provides a 
potential lower bound to F . In particular, when the estimate for Fmin obtained from Eq. ( p7|) 
exceeds Gmax i it is quite possible that pair creation never occurs at all, since opacity for photons 
of energy £max occurs through interactions with photons at even higher energies, for which there 
is no observational evidence. Cases with values of Fmin/emax greater than unity generally arise 
only for small At and at cosmological distances (see Table |3|: GRB910601 is a perfect example). 
In such instances, these cutoff considerations become quite relevant, and it becomes necessary to 
take Emax as the estimate for the minimum Lorentz factor F . Since values of Emax between 1 
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and 10 MeV in the collection of EGRET bursts are a marker of fainter or steeper spectrum bursts, 
i.e. probably reflecting the observational limitations of EGRET, it is quite realistic to develop bulk 
motion estimates based Eq. (p7[), in the belief that many (if not most) bursts emit at energies much 
higher than 100 MeV. 



3.2. Generalized Geometries 

Generalizing from Eq. ( |47| ) to the expansion geometries considered in this paper, a pair produc- 
tion transparency condition is obtained by setting the optical depth that is obtained from Eqs. ( ^7[ ) 
and ( [T7| ) equal to unity, i.e. effectively reading off abscissa values for chosen ordinates from curves 
like those in Fig. ^. Since the expressions for the optical depth and flux from the expansion involve 
integrals with integrands that depend on the bulk Lorentz factor, the solutions for Fmin roots to 
r^^ = 1 must be solved iteratively: we adopt a bisection technique. The minimum bulk Lorentz 
factors Fmin for EGRET burst sources that result from our pair production transparency calcu- 
lations in Eqs. (|37|) and (^) are depicted in Fig. ^ as functions of the fractional shell thickness 
AR/Ro . This illustration limits shell thicknesses to regimes where AR/Rq < 1 , since Rq is tied 
to the time variability via Rq = TcAt , and this coupling becomes inappropriate for AR ^ Rq 
regimes, where AR ~ cAt is a more apt choice. Again the source parameters from Table ^ are 
used, and results are presented for large expansion half-angles, 0b = 90° , variability timescales of 
At = 1 ms, and source distances of 100 kpc and 1 Gpc that represent galactic halo and cosmological 
burst scenarios, respectively. Cosmological redshift modifications, which depend on the choice of 
cosmology, are neglected for simplicity. All of the twelve curves in the four panels exhibit similar 
behaviour, with F^in being independent of 5 = AR/Rq when 5^1-/3, the so-called thin-shell 
limit, and F^in declining roughly as for 5 ^ 1 — (5 when the filled-sphere regime is realized. 



These dependences on AR/Rq appear explicitly in the asymptotic forms in Eqs. (41) and (42). 
Domains where AR/Rq > 1 (i.e. for a true filled sphere) are not depicted since they are unlikely to 
be encountered in gamma-ray bursts. The values of F^in obtained when AR/Rq ~ 1 — /3 (i.e. the 
transition regions) are comparable to those listed in Table |3| for all bursts; this reflects the broad 
applicability of the constraint in Eq. (^) that is the hallmark of the so-called "blob" calculation. 

The curves in Fig. ^ generally concur with the global trends of increasing F^in for higher 
Emax and/or declining spectral index a. Yet the cross-over of two curves in the lower left-hand 
panel exemplifies how the expansion geometry can complicate trends, and generate non-monotonic 
behaviour in Smax or a . For comparison, Fig. Q reproduces the bottom right panel of Fig. ^, but 
for a variability timescales of At = 1 sec. The curves resemble those of Fig. ^. However, at these 
cosmological distances, F^™ is reduced from the corresponding values in Fig. ^ solely because the 
large At dilutes the density of internal photons inferred for the source. Again a cross-over of curves 
arises, indicating that non-monotonicity of F^in in Emax and a does not belong exclusively to 
galactic halo scenarios. Note that the uncertainties in the observational quantities that are listed 
in Table are as large as around 10%. These lead to uncertainties in the Fmin determinations for 
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both these figures of the order of 20%, which are largely masked by the ranges of Tmin produced 
by varying model parameters. Hence further consideration of experimental uncertainties is omitted 
from this paper. Note also that the results we have presented have assumed that the test photons 
in the source start from the rear of the expansion. This maximizes the computed optical depths, 
implying that it is possible to lower our estimates for Fmin somewhat. Permitting test photons 
to be emitted throughout the source will perhaps lower the mean optical depth by a factor of 2 
or so, leading to a reduction of r^m of around 15%. Hence the detailed consideration of the 
distribution of test photons is neglected in this paper as this will only have a minor impact on the 
Tmin obtained. 

The results presented in Figs. ^ and § clearly define the behaviour of our Fmin solutions 
for large opening angles, but obviously do not address the entire available phase space for model 
parameters. Hence in Fig. |5| we depict the variations of Fmin with expansion opening angle 0b 
for fixed AR/Rq. These curves unequivocally demonstrate how insensitive T^in is to 0b when 
0B ^ 3AR/Ro , i.e. for a range from modestly small 0b up to 90° . Note that such a substantial 
range is realized for the adopted values of AR/Rq , chosen to correspond to the transition between 
thin-shell and thick-shell regimes. This insensitivity, whose important implications for cosmological 
models are discussed at the end of this subsection, is a principal conclusion of this paper. The values 
of Fmin ^-t the intercept with the right hand side axis of each panel in Fig. ^ are just those obtained 
by vertically slicing the plots in the bottom two panels of Fig. ^ at the appropriate value of AR/Rq . 
For 0B ^ AR/Rq, power-law asymptotic behaviour of Fmin is observed in the Figure, with the 
dependence being easily deduced (for 1 — /? <C AR/Rq ) from Eq. (|43|): Fmin oc 
is a very weak dependence (for a in the range 2-3 typical of EGRET bursts) on 0b , generated 
by the strong variation of the optical depth with F . 

It is quite instructive to augment these plots by summarizing the behaviour of the optical 
depth results with phase space diagrams using well-chosen variables. This can be achieved in an 
enlightening manner for both observational and theoretical (i.e. model) phase space parameters via 
contour plots, i.e. exhibiting curves of constant F ( = Fmin ) that satisfy the criterion r^^ = 1 for 
the pair production optical depth. First of all, we focus on the space of observational parameters 
given by the maximum energy observed, Emax , and the EGRET spectral index a , both being 
listed in Table |2| for EGRET bursts. Fixing the observed flux and the source time variability at the 
"canonical" values of /(I MeV) = 3cm~^ sec" -1 MeV"^ and At = I ms, respectively, the resulting 
contour plot is shown in Fig. These contours, which represent lower boundaries to regions of 
opacity (i.e. r^^ > 1 )> display a number of trends that are hallmarks of the optical depth properties 
of the relativistically expanding radiation gas. 

First, Emax is an increasing function of a . This arises because, for the particular values of 
/(I MeV) , At and F chosen, solutions with Emax < F are always realized. For such solutions, 
the test photons at energy Emax interact with photons near pair production threshold in the CM 
frame, i.e. with photons of energy around Emax/F^ ( < 1 MeV) in the observer's frame. Since the 
optical depth is held constant (i.e. unity), and the flux is pinned at 1 MeV, thereby providing a 
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"pivot point" in the spectrum, increasing a then raises the number of interacting photons (with 
energies below 1 MeV) so that Emax must correspondingly be increased to compensate. This trend 
solely finds its origin in the realization of an e^Ax < T branch of solutions. As q ^ oo , the 
energy of the interacting photons must approach the pivot point, i.e. 1 MeV. Hence the Smax 
curves asymptotically approach ~ as a becomes very large, behaviour that is conspicuous 
in the galactic halo cases in Fig. ^. In the particular examples shown, Smax drops off rapidly 
as a approaches unity, a singularity of these curves [e.g. see Eq. (^)]. The monotonic increase 
of £max with a can be inverted to yield a declining tendency if an Smax > T solution branch 
can be encountered, so that interacting photons are always above the pivot point at 1 MeV. This 
occurs when the product /(I MeV) d^/At is small enough, i.e. forcing Smax higher for given F 
and Q [e.g. see Eq. (^7|)]. Then the contours would rise rapidly as a approached unity, but 
still asymptote to roughly a F^ dependence for a ^ cxd . This situation is more likely to arise for 
At = 1 sec variability timescales. Note that the larger values of /(I MeV)(i^/At in the cosmological 
cases in Fig. ^ yield stronger dependences of Emax on a , primarily because these cases generally 
have larger "lever arms" for the interacting photons (at energy ~ ^max/E^ ) around the pivot point 
at 1 MeV. 

The dependence of Emax on F is very closely given by the blob calculation in Eq. (^7|), 
namely Smax oc r-{i+2a)/(a-i) _ por example, doubling F at a = 2 in either of the galactic halo 
or cosmological cases yields an increase of around 29.4 , which is close to the "blob" estimate of 
32 . Similarly, doubling F at a = 3 in either of the galactic halo or cosmological cases yields 
an increase of around 11.2, which is very close to the "blob" estimate of 11.3. For fixed a, 
these amplification ratios are independent of the source distance d, since d just forms part of 
the proportionality constant for the relationship between Smax and F . Another trend that is 
apparent in Fig. ^ is that Smax declines with decreasing b = AR/Rq . The variation of Smax with 
6 is depicted using the light solid curves for the galactic halo scenario, with the 5 = 0.3(1 — (3) 
and 6 = 3(1 — /?) cases visually defining a band around the 6 = 1 — P case (the behaviour for 
cosmological source distances is similar). When 6 is reduced, the average density of photons within 
the source increases, pushing the optical depth up. Hence, to compensate, Emax must also decline 
with S , so defining the observed trend; in Figs. |3| and ^, this effect forces F^in to increase when 
Emax is held constant. From those figures, it is evident that when (5 <C 1 — /3 , the thin-shell limit 
produces insensitivity of the optical depth to (5 , a feature that is also apparent in Fig. ^, for which 
the 5 = 0.3(1 — P) curves are more proximate to the (5 = 1-/3 ones than are the 5 = 3(1 — /3) 
cases. Finally, we note that the value of /(I MeV) for each of the EGRET bursts that are depicted 
as points in Fig. ^ differs from the chosen canonical value. Hence for those sources with higher 
/(I MeV) (GRB 910503 and ORB 910814), the exhibited curves should be slid down somewhat to 
visualize the situation (i.e. infer bulk Lorentz factors) appropriate for these bursts. Likewise, for 
the remaining EGRET bursts, the curves should be moved upwards to deduce F values that are 
consistent with the results depicted in Figs. ^ and |^. 
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The theoretical phase space contour plot is presented in Fig. |^ which exhibits curves of constant 
r that satisfy the criterion r^^ = 1 , in the space defined by the opening angle 0b and the fractional 
shell thickness 6 = AR/Rq . We focus on regimes where AR/Rq < 1 since, as mentioned above, 
we have tied Rq to the time variability via Rq = TcAt , a coupling that becomes inappropriate 
for AR ^ Rq regimes, where AR ~ cAt is a more apt choice. The large range of ©b is chosen 
intentionally to present the information relating to the reduction of expansion opening angles that is 
omitted from earlier considerations, such as in Figs. Hjand^ For all contours, the maximum energy 
and spectral index were set at Smax = 100 MeV and a = 2.21, respectively, and the variability 
timescale At and the flux /(I MeV) at 1 MeV are representative of EGRET- like sources. These 
parameters were tuned somewhat to obtain a maximum of informational content in the figure, so 
that other choices of parameters can lead to some variations in contour shape. Remembering the 
general trend of a reduction in r-y-y with increasing 0b , the contours in Fig. |^ clearly represent 
upper boundaries to regions of opacity (i.e. r-y-y > 1 ). 

A number of prominent features appear in this phase space plot. Foremost among these are 
the vertical portions of the F = 10 and F = 250 contours, present when the opening angles are 
significant. These define regimes where the optical depth is independent of ©b , and both thin-shell 
and thick-shell regimes that are well-described by the asymptotic formulae in Eqs. (^) and ( ^2|) 
can be realized. It is precisely this upper region of the 0b - S diagram that is probed in solutions 
depicted in Figs. |^ and The virtual independence of the optical depth to 6 observed in those 
solutions for F^m manifests itself in Fig. |^ as an extreme sensitivity of the position of any vertical 
sections of the contours to the choice of F (or Emax or a). This sensitivity therefore produces a 
low density of contours in the upper left hand portion of Fig. |^, so that contours possessing vertical 
sections occupy a minority of cases if broad ranges of ©b are considered. Note also that low values 
of r are obtained only in the upper right of the figure. However, since values of AR/Rq greater 
than unity are unrealistic, it becomes clear that for this choice of Emax and a , only values of 
F ^ 5 are attained. This signifies the general property of these calculations that relativistic bulk 
motions are always inferred unless Emax is not too much greater than 1 MeV. 

The narrow beam (i.e. ©b <^ 0.01 ) portion of the parameter space exhibits distinctive power- 
law dependences, with ©B rising as \/J when 6 is very small, and declining as 1/6 when 6 exceeds 
1 — (5 . This asymptotic behaviour can be deduced with the aid of Eqs. (H) and (||). Table | 
identifies four parameter regimes, three of which are relevant to narrow beam considerations, the 
other being the domain of large ©b discussed in the previous paragraph. Eq. ( ^4|) is pertinent to 
the lower left hand portion of Fig. 0; contours of unit optical depth thereby defining the dependence 
©B oc (5i/2/r<^+i/2 (^for Rq = TRy). In this limit, doubling F decreases ©b by around a factor 
of 6.54 for a = 2.21 , regardless of the assumed distance to the source. The lower right of the 
figure is described by the thick shell (i.e. 5 ^> 1 — /?) limit of Eq. (|4^), which yields contours 
with ©B oc (5~^F~(^'^"'"^) ; doubling F in this limit reduces ©b by a factor of 343 for a = 2.21 , 
behaviour that is borne out in Fig. |^. The third asymptotic domain is defined by the thin-shell 
limit of Eq. (|^), yielding contours approaching a limit with ©b cx r^(^"+^) , independent of 5 . 
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This domain is almost attained at the broad peaks of the contours that remain always in the narrow 
beam (i.e. lower) portion of phase space, particularly for the F = 40 and T = 1000 cases in the 
figure. These three limiting forms can be written explicitly as (for F ^ 1 and Rq = TR^ ): 

( ^2 njf \ ^ 1/2 {a-l)/2 

These asymptotic formulae are depicted as thin, light, dotted lines in Fig. for the F = 40 galactic 
halo case, clearly indicating how the contours closely approach these in the appropriate ranges of 
6 . The lowest F examples in Fig. for each of the galactic halo and cosmological scenarios do not 
realize thin shell portions of phase space without assuming significant opening angles 0b . This 
feature marks the general property that low F curves occupy the upper right corner of the 0b -5 
diagram, a domain where the thick shell parts of the solutions in Figs. |3| and ^ are appropriate. 
Other trends relevant to the 0b diagram include a general reduction of 0b with increases in 
a or decreases in Emax • 

This concludes the survey of observational and theoretical parameter space. In view of the 
extensive presentation of results in this section, it is important to highlight the implication of this 
work that is most salient for gamma-ray bursts. The principal conclusion of our work (expounded 
in brief in Baring and Harding 1996, with a preliminary version given in Harding and Baring 1994), 
is that Fmin is quite insensitive to the choice of 0b when 0b ^ 1/Fmin , behaviour that can be 
inferred from Figs. |2[ This result arises because causality restricts the available phase space for 
pair production interactions more effectively than does the expansion opening angle 0b , when 
0B ^ 1/Fmin ; it has profound repercussions for gamma-ray burst models. For such source models, 
the principal advantage (e.g. Krolik and Pier 1991) of restricting 0b to small values like 1/F 
is a lower (solid angle-reduced) luminosity at the source for a given observed flux. However, the 
number of non-repeating sources must then be a factor 0~^ ~ F^ higher in order to account for the 
observed burst rate. In the case of cosmological GRBs, this factor could be as high as 10® for the 
values of Fmin determined here, which is unacceptably large for many models, particularly those 
that involve neutron star-neutron star or neutron star-black hole mergers (Paczyhski 1986; Eichler 
et al. 1989; Narayan, Piran and Shemi 1991; Meszaros & Rees 1992), failed Type lb supernovae 
(Woosley 1993) and rapid spin-down of high- field millisecond pulsars (Usov 1992). This defined the 
commonly-perceived "number problem" for cosmological bursts. Clearly, in view of the results of our 
analysis, this problem is a non-issue, since imposition of small opening angles 0b in order to satisfy 
pair production transparency in EGRET bursts is not necessary. Hence, causality restrictions to the 
optical depth differ so little between 0b = 90° and 0b ~ 1/F cases that burst population statistical 
requirements can comfortably be satisfied without resorting to beamed expansion geometries. Of 
course, opening up the expansion angle then amplifies the energetics requirements accordingly, so 
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that model development must still meet the needs of acceptable bursting rates and energy budgets. 
In such considerations, it is evident from the work presented here that pair production constraints 
will play only a secondary role in determining such model requirements, becoming involved purely 
through the evaluation of permissible Lorentz factors for bulk motion in gamma-ray bursts. 

3.3. Discussion 

The choice of coupling the scale of the expansion to the variability timescale via transverse 
dimensions is subjective, though it is widely adopted in applications of bulk relativistic motion in 
astrophysics. Other choices are possible, such as using burst or subpulse durations (e.g. Fenimore, 
Madras and Nayakshin 1996) and/or relating these to longitudinal dimensions in the source. It is 
fitting to outline the reasons for adhering to our preference. Suppose that timescales larger than the 
variability time At , for example the burst duration , are used as the observational diagnostic 
of the source size. For uniform expansions, as we have assumed, the geometrical appearance of the 
"look-back" volume forces the time profile to maintain a well-defined shape ( [t/Td\~°'~'^ for spectral 
index a = 2) that mimics the so-named FRED (fast rise, exponential decay) profile (Fenimore, 
Madras and Nayakshin 1996). This profile necessarily has a width of the order of under these 
assumptions, so that its smooth, decaying shape is inconsistent with the vast majority of burst time 
histories. Temporal consistency can therefore be attained only if the appropriate observational 
timescale is of the order of the variability time 6t , so that the burst comprises a multitude of 
shells, or perhaps if a single shell is "patchy" in the transverse dimension. The latter possibility 
still produces time profiles that do not match many burst histories, so that one is compelled to 
adopt the many-shell proposition, perhaps produced by a central engine, as advocated by Fenimore, 
Madras and Nayakshin (1996). The timescale that is then appropriate is At , precisely our choice, 
though the value of this depends on whether BATSE or EGRET variabilities are used (as discussed 
in Section 3.1 above). 

The issue of whether the variability should be tied to transverse or longitudinal source dimen- 
sions remains to be addressed. Since At is always close to the threshold of temporal resolution 
of any of the CGRO instruments BATSE, Comptel and EGRET, it is appropriate to assume that 
measured variability is actually an upper bound to the source variability. If we opt to relate this to 
the direction transverse to the line of sight to the observer, then the inequality Rq/^ ^ cAt follows, 
and it is customary to take the equality to specify Rq . By the same token, if dimensions along the 
line of sight are preferred, this inequality is replaced by Rq/T'^ ^ cAt . A consistent description of 
the expansion can only be obtained when both of these inequalities are satisfied, which obviously 
occurs when the more-constraining Rq/T ^ cAt is adopted. This motivates our choice of coupling 
the variability to the transverse dimension; tying it to the line-of-sight direction is insufficiently 
restrictive. Notwithstanding, the difference between these two choices is merely one Lorentz factor 
in the optical depth [compared with around five or six imposed by the spectrum: see Eq. (^3[)], 
to which the estimates of Fmin are quite insensitive: opting for Rq/T^ <^ cAt reduces Fmin by 
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factors of the order of two or less. Other possibihties for choosing the scale of the expansion exist, 
such as AR = cAt , and these are discussed at length in the temporal analysis of Fenimore, Madras 
and Nayakshin (1996). Note that fixing AR = cAt with either Rq ~ FcAt or Rq ~ T^cAt yields 
AR/Rq <C 1 , comfortably in the phase space covered by Fig. |^. The essential point that should 
be emphasized is that these subjective alternatives probe details of the expansion microstructure 
that are beyond the purpose of this analysis, and are largely peripheral to it, primarily because 
of the relative insensitivity of the F^in estimates to these choices. The principal conclusions of 
this paper, including the insensitivity of the optical depth to the expansion opening angle 0b , are 
guaranteed regardless of such variations on our assumptions. 

One question that naturally arises when obtaining estimates for the bulk Lorentz factors via 
pair production constraints is why the values of 100-1000 obtained here for cosmological bursts are 
of the same order as those obtained from fireball expansions (e.g. Paczyhski 1986; Shemi & Piran 
1990; Rees &: Meszaros 1992) of enormous initial optical depths. This similarity is no coincidence. 
The bulk motions attained by the adiabatic expansion phase of pure electron-positron (pair) fireballs 
yield Lorentz factors F that saturate at some value corresponding more or less to the "freeze-out" 
of pair production, i.e. the epoch of free expansion is approximately marked by the onset of pair 
production transparency. For cosmological bursts, where the luminosities can be of the order of 
L ~ 10^^ — lO^'^ erg/sec and the energy deposition can be larger than 10^^ ergs, the optical depth is 
roughly La^ / {RmeC^) /T^ (e.g., for E~'^ spectra) and can be 10^°/iZio or larger for F = 1 , where 
RiQ is the size of the fireball at the end of the epoch of opacity in units of 10^'^ cm. This optical 
depth can be reduced to unity by relativistic beaming with F in the range of 100-1000 when 
i?io ~ 1 . Since the Lorentz factor attained during the fireball "acceleration phase" scales roughly 
as its radius (e.g. Paczyhski 1986; Piran, Shemi and Narayan 1993), then it follows that values of 
a freeze-out radius of Riq ~ 1 would correspond to F ~ 100-1000 for fireballs initiated in regions 
of diameter 10^ - 10^ cm. These lO^'^ cm scalelengths for the onset of expansion transparency are 
comparable to those used for Rq = TcAt in this paper, thereby explaining the similarity of our 
estimates for Fm™ to the Lorentz factors of fireball- initiated relativistic expansions. 

The results we have presented focus on the energy range appropriate to EGRET detections 
of gamma-ray bursts. There are now ongoing programs for searches for bursts at TeV energies, 
specifically the target of opportunity monitoring of BATSE localization error boxes by the Whipple 
air Cerenkov experiment, using the rapid response that is facilitated by the BACODINE alert 
network (Barthelmy et al. 1995). While these efforts have failed to provide any positive TeV 
detections so far, probably due to the fact that Whipple's sensitivity threshold still inhibits any 
possibility of detection for all but the very brightest of bursts (see Connaughton et al. 1995 for a 
discussion of the current Whipple sensitivity), the prospect of large field of view monitoring of the 
sky by the air Cerenkov water tank detector MILAGRO (e.g. Yodh 1996) in the very near future, 
promotes the extension of our bulk motion estimates to the TeV energy range. Such considerations 
also anticipate future space missions like GLAST, which will span the 10 MeV-200 GeV range. 
Obviously, increasing Smax to TeV-type energies would tend to push estimates of the bulk Lorentz 
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factor up, in order to suppress pair creation. To explore the implications of TeV-emitting bursts 
sources for estimates of r^m > we computed the infinite power-law "blob" calculation solutions to 
Eq. (^7[), and depicted them in Fig. ^, as a function of the spectral index a/j . It is sufficient to 
focus on this simplest of cases, noting that the complicating effects of expansion geometry mirror 
those considered at lower maximum energies. 

In the figure, the cosmological cases exhibited the expected trend of a dramatic increase in 
Fmin for flatter spectra, a consequence of the enhanced supply of interacting photons at ~ F^/smax 
(generally well above 511 keV) for lower values of . For typical EGRET source spectral indices, 
in the range 2-3, F^in is indeed an increasing function of Emax • However, for a/i < 1 , this 
behaviour is reversed, with Fmin declining with Emax , because the optical depth in Eq. (p7|) is 
then a decreasing function of Smax • Note that in the figure, the galactic halo case ( d = 100 kpc) 
displays a comparative insensitivity of Fmin to ah ■ This insensitivity arises because the energies 
( ~ F^/e^Ax ) of the photons interacting with those at Emax are generally relatively near the pivotal 
energy of 511 KeV, where the source flux is pinned. In Fig. ^ the source flux / at 511 keV is typical 
of BATSE burst detections; for this flux, the MILAGRO experiment will be sensitive to bursts with 
^ 2.6, a dividing line that is marked in the figure. Hence, only bright, flat spectrum sources 
like GRB 910503 provide good candidates for potential detections at TeV energies by MILAGRO, 
and for that matter Whipple. GLAST will have the capability of spanning the EGRET and sub- 
TeV energy ranges. This extension of F^m estimates to Smax values in the TeV range of course 
assumes that bursts intrinsically emit at such high energies. Note that for the cosmological case, 
such intrinsic burst emission is also subject to pair creation in collisions with photons supplied by 
radiation fields external to the burst, for example the infrared background (Stecker and De Jager 
1996, Mannheim, Hartmann and Funk 1996). This type of attenuation is strongly dependent on 
the redshift of the source, so that sources at 1 Gpc would be strongly attenuated at 1 TeV, while 
those at 100 Mpc would be transparent. Such external attenuation considerations are beyond the 
scope of this paper, and are examined in detail by Mannheim, Hartmann and Funk (1996). 

The analysis of this paper has made the expedient assumption that the burst spectra are infinite 
power-laws, with the spectrum matching that observed by EGRET. This, of course assumes that the 
turnovers seen at sub-MeV energies in BATSE data are immaterial to pair production calculations. 
The relevance of spectral curvature below the EGRET range to attenuation studies depends upon 
the energies of the photons that interact with the test photons at et = £max , which generally are 
around F^/smax in the stationary observer's reference frame for infinite power-laws, as discussed 
just above. Clearly, low energy spectral paucity (relative to EGRET-range power-laws) at around 
r^/^MAx limits contributions near threshold so that phase space near the pair production threshold 
is inaccessible, and the process is pushed into the Klein-Nishina regime: the optical depth drops 
accordingly. Hence MeV and sub-MeV spectral curvature plays a role in opacity determinations 
when F^/e^Ax ^ 1 1 a situation that manifests itself at galactic halo (or even galactic disk) distances 
for most of the EGRET bursts, as is evident from an inspection of Tables § and ^. Consequently, 
regimes of transparency may become possible at high energies in galactic bursts, so that the range 
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of opacity becomes finite. These issues are discussed in detail in Baring and Harding (1997), where 
it is demonstrated that reaUstic broad-band GRB spectra may yield broad absorption troughs 
(and also other spectral forms such as shelfs) in the 1 GeV - 1 TeV range of bursts spectra in 
galactic halo sources, distinctive spectral features whose existence or otherwise could be probed 
by experiments such as GLAST, Whipple and MILAGRO. Note however, that for cosmological 
sources, the inferred Fs are so high that F^/emax always exceeds unity and the infinite power-law 
calculations presented here are always quite appropriate. As a result, the optical depth is then 
always monotonically increasing in Emax ) so that only simple (exponential or power-law) spectral 
turnovers are possible. Hence, Baring and Harding (1997) postulated that such distinguishable 
spectral structure may provide a means of discriminating between a galactic or cosmological origin 
for gamma-ray bursts, an enticing prospect for the high energy astrophysics community. 

4. CONCLUSION 

In this paper, we have presented our extension of pair production transparency calculations 
in relativistically expanding gamma-ray burst sources to quite general geometries, including shells 
of finite thickness and arbitrary opening angle. This work includes an extensive analytic reduction 
of the optical depth from a quintuple integral to a single integral in the special, but quite broadly 
applicable, case of observing photons only along the axis of the expansion. Such a reduction 
is extremely expedient for opacity and transparency considerations, providing the reliable and 
numerically-amenable analytic expressions in Eqs. (|l7|) and ( ^7[) that completely describe the pair 
production optical depth. We determine that the minimum bulk Lorentz factor Tmin for the 
EGRET sources to be optically thin up to the maximum energies observed, i.e. display no spectral 
attenuation, is only moderately dependent on the shell thickness and virtually independent of 
its opening solid angle if 0b ^ I/Fmin- This insensitivity to ©b, which is a consequence of the 
strong impact that causality has on the available interaction phase space, relieves the commonly- 
perceived number problem for non-repeating sources at cosmological distances: it is not necessary 
to invoke small ©b to effect photon escape. This negation of the number problem for a wide range 
of expansion geometries is the principal conclusion of this paper, and is an important result for 
specific cosmological burst models. The values of F obtained, typically of the order of 10-30 for 
halo bursts and ^ 100 for sources of cosmological origin, depend only moderately on the choice of 
GRB timescale used to determine the expansion size. Our new limits on required expansion velocity 
for given source geometries will significantly aid the placing of realistic constraints on gamma-ray 
burst source models. 

We thank Brenda Dingus and Jennifer Catelli for many discussions about EGRET burst data, 
and Ed Fenimore and Jim Ryan for numerous conversations concerning gamma-ray bursts. This 
work was funded, in part, by the Compton Gamma-Ray Observatory Guest Investigator Program. 
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A. APPENDIX 

PROPERTIES OF THE HYPERGEOMETRIC FUNCTION Ga{z) 



In the integrand of the expression for the pair production optical depth, Eq. (|37D, appears the 
hypergeometric function 

Ga{z) ^ dq = + 1. 1; 2a + 2; z) . (Al) 

This can be represented by the alternative hypergeometric form 

= TT2^T^^(« + ^' ^' ^" + ^' ' ^^'^^ 

using the transformation formula 9.131.1 of Gradshteyn and Ryzhik (1980). 

The argument of Qa{z) is either A or crA , which can be derived from Eq. (^6|). It is clear that 
for s > 1 , A attains values within a range — oo < A < Amax , with A^ax = + /?) (note that 
dX/ds > 0). Clearly the condition C < 1 renders a less than unity, so that crA is also bounded by 
the range (— oo, /3/(l + /?)] . It follows that there are two natural ways to evaluate Qa{z) for the 
purposes of this paper, both using the series expansion (9.100 of Gradshteyn and Ryzhik, 1980): 

Fia + 1 V2a + 2- z) - 1 + " + ^ z + (« + !)(« + 2) ^2 .^os 
l^{a + I, I, la + I, z) - -L + 2a + 2 ^ + (2a + 2)(2a + 3) ^ + ' ' ' ' ^^"^^ 

When \z\ < 1 , this series can be used directly with convergence as rapidly as the geometric series 
^^(z/2)" . When — oo < z < — 1 , the alternative form for Qaiz) in Eq. ( |A2D can be used, where 
1/2 < z/{z — 1) < 1; the series in Eq. ( |A3D with the substitution z z/{z — 1) then also 
converges like the geometric series J2niz/'^[z — 1])" j i-e. with the same rapidity. With this scheme, 
computation of Gaiz) to high accuracy is quick. 

Note that as s ^ 1, A — > — oo. Using the representation in Eq. ( [A2| ) and also Eq. (36), it 
is clear that in this limit (1 + 2a)^a(A) approaches F{a + 1, 1; 2a + 2; 1)/(1 — A) , and therefore 
becomes approximately proportional to s — 1 . It follows that the integrand in Eq. ( |37|) is finite as 

Also of use in this paper, specifically in the determination of the optical depth in the limit of 
filled spherical expansions, is the integral identity 

\zzgUz) = ^ + ^{^{2a)-^ia)-l} , ^(x) = ^{log,r(x)} , (A4) 

for use in /5 « 1 situations. This can be established using the integral representation of Ga in 
Eq. ( [AID , then reversing the order of integration and performing the z -integration analytically. 
An integration by parts then enables the use of identity 3.231.5 of Gradshteyn and Ryzhik (1980), 
and the result ensues. The finite rational series for ij;{x + n) — ^(x) was also used in manipulating 
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Eq. ([A4|), and for integer a it can be used to obtain rational values for these integrals. For /3 <C 1 
cases of initially filled spherical expansions, the integral 



dz 



Qa{z) 



■k/2 



dO 9 cos^" I 



1 



l + 2a 



1 + (vr - 3) 



2(a + l) 



9/8 



(A5) 



is needed. The identity is established by using the transformation 9.134.1 in Gradshteyn and 
Ryzhik (1980) for the hypergeometric function in Eq. ( [All) , changing to an integration variable 
of [z/{2 — z)]^ , using the integral representation for general hypergeometric functions in 9.111 of 
Gradshteyn and Ryzhik (1980), and then reversing the order of integration. The integral, when 
multiplied by 1 + 2a , is only weakly dependent on a , and the simple approximation obtained in 



Eq. (A5) is accurate to better than 1% for < a < 1/2 and better than 0.1% for 1/2 < a < 10 . 
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TABLE 1 

Optical Depth Dependence on 5 and 0b for Narrow Beams 



Regimes of S 
for 0B < 1 - /? 


Eq. (17) 


Eq. (37) 


"^77 


0| <0B < 
0| <0B <^ < 1-/? 
0| <^ <0B < 1-/? 
(J < 0| < 0B < 1 - /? 


^0| 

^0| 
J2 


^2 0B 
^ 0B 

5^ 

J2 


1 

^0B 
1 
0B 

el 

constant 



Note. — The proportionalities, of the optical depth t^^ and 
the two factors T j H and r^^j H contributing to it, to the expan- 
sion parameters J = I^^R/Rq and ©b, specifically in the limit of 
narrow beam expansions. Remember that in the text, ^ = cos ©b 
so that 1 — C ~ 0|/2 , and various regimes of shell thickness or 
thinness are described in the left-hand column. 
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TABLE 2 

Parameters for Energetic Gamma-Ray Bursts 





Detecting 






/(I MeV) 


GRB 


Instrument 


(MeV) 


a 


(cm-^sec-^MeV"^) 


910503 


EGRET 


170 


2.2 ± 0.2 


8.71 ± 0.49 


910601 


EGRET 


3.5 


3.7±0.2 


0.98 ± 0.08 


910601 


COMPTEL 


5.0 


2.8 ± 0.2 


0.50 ± 0.10 


910814 


EGRET 


60 


2.8 ± 0.2 


13.5 ± 0.85 


930131 


EGRET 


1000 


2.0 ± 0.3 


1.95 ± 0.26 


940217 


EGRET 


3380^ 


2.5± 0.1 


0.36 ± 0.03 


950425 


EGRET 


120 


1.93 ± 0.04 


1.62 ± 0.09 



Note. — The maximum energy Smax of detection, the high energy spectral 
index a , and the (measured or extrapolated) source photon flux at 1 MeV for 
six bursts observed by EGRET. The EGRET detections are from Schneid et al. 
(1992: GRB 910503), Kwok et al. (1993: GRBs 910601 and 910814), Sommer 
et al. (1994: GRB 930131), Hurley et al. (1994: GRB 940217) and CateUi et 
al. (1996 GRB 950425). The COMPTEL measurement of GRB 910601 (see 
Hanlon, et al. 1994) is included because it was seen at a higher energy than the 
EGRET detection of this source, and is used in the estimates for lower bounds 
to r in Table 3. -"-The famous 18 GeV photon from GRB 940217 was omitted 
because it was not contemporaneous with any emission below 100 MeV. 
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TABLE 3 

Lower Bounds for T in Energetic Gamma-Ray Bursts 



GRB 


At 


1 kpc 


100 kpc 


P 

■L MIN 

100 Mpc 


1 Gpc 


JVlciximiini distciiicc 
for isotropic emission 


910503 (E) 


1 ms 


2.9 


16 


199 


460 


49 pc 


910601 (C) 


1 ms 


1.0 


3.8 


31 


61 


1160 pc 


910814 (E) 


1 ms 


3.1 


12 


101 


203 


24 pc 


930131 (E) 


1 ms 


2.8 


17 


276 


694 


79 pc 


940217 (E) 


1 ms 


4.2 


19 


195 


419 


13 pc 


950425 (E) 


1 ms 


1.6 


11 


188 


486 


297 pc 


910503 (E) 


1 s 


1.0 


4.5 


56 


130 


1.6 kpc 


910601 (C) 


1 s 


1.0 


1.3 


11 


22 


37 kpc 


910814 (E) 


1 s 


1.1 


4.4 


35 


71 


0.76 kpc 


930131 (E) 


1 s 


1.0 


4.4 


69 


174 


2.5 kpc 


940217 (E) 


1 s 


1.3 


6.2 


62 


133 


0.43 kpc 


950425 (E) 


1 s 


1.0 


2.6 


45 


117 


9.4kpc 



Note. — The estimated minimum bulk Lorentz factor Fmin from Eq. (47), and the maximum 
source distance that permits isotropic emission, for the six GRBs detected by EGRET, assuming 
that two-photon pair production does not influence the spectrum below the maximum detected 
energy £max • The observed spectral index a below ^max , as well as the source power-law flux /(I 
MeV) at 1 MeV are given in Table 2 (E labels use EGRET data and C indicates use of COMPTEL 
data) ; the EGRET power-law fits were extrapolated down to 1 MeV where necessary to define the 
normalization measure /(I MeV). At 1 kpc, Fmin = 1 entries represent Fmin < 1 solutions to (i.e. 
the breakdown of) Eq. (47). Estimates for each burst are given on the two timescales At typical 
of source variability and subpulse duration, as seen by BATSE. Note that source redshifts are set 
to z = in order to avoid a choice of cosmology. 
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Fig. 1. — (a) A depiction of the source geometry for the relativistic expansions of bulk Lorentz 
factor r considered in this paper. The emission region is at aU times a conical sector of a spherical 
shell with half-angle 0b • At some time t the inner radius of the shell is given by R , while 
the thickness of the shell is always AR . The test or observable photon is marked by T, and it 
can interact via 77 — > e^e~ with other (so-called interacting) photons at typical position I. The 
angles their position vectors (from the origin O of the expansion) make with the axis OZ of the 
expansion are O* and 0^ , respectively; the angle subtended by these vectors at the origin is Qu ■ 
(b) A depiction of the spatial and angular variables, as defined in the text, that are relevant to the 
analysis of pair production in expanding burst sources. OZ represents the axis of the expansion 
and generally is not co-planar with the plane (OTI) formed by the positions of the test (T) and 
interacting (I) photons and the origin (O) of the expansion. The two photon momentum vectors 
(heavy lines) have an angle 9ti between them and define a third plane that is generally not co- 
planar with OTI. The "angles of non-radiality" of the test and interacting photons are 6t and 6i , 
respectively, and Ou defines the angle subtended at the origin by the positions of the two photons. 
The labelled distances are rt = \OT\ , = \OI\ and ru = \TI\ . 
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Fig. 2. — The optical depth T^^{et) of test photons to pair production, scaled as in Eq. (^6|), 
but for redshift z = , as a function of bulk Lorentz factor T for different opening angles 0b 
of the expansion, as labelled, and for shell fractional thicknesses (a) AR/Rq = 10~^ and (b) 
AR/Rq = 10~^ . The curves are obtained for a typical EGRET source spectral index of a = 2.5 
and for Rq = TcAt (discussed in the text). The scaling distance is chosen to illustrate typical 
optical depths of cosmological bursts, however the curves are easily translated down by around 8 
orders of magnitude if a galactic halo hypothesis is preferred. 
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Fig. 3. — The minimum bulk Lorentz factor Tmin for six EGRET GRBs, as obtained from the 
pair production condition T^^{et) = 1 in Eq. (37). Here et is the maximum energy Emax detected 
by EGRET; values of £max cind other observational parameters are listed in Table Q (Comptel 
data are used for GRB 910601, as in Table ^). The top two panels consider the first three EGRET 
bursts, while the bottom three are for the most significant of more recent events. Results are shown 
for two different source distances, d = 100 kpc (left-hand panels) and d = 1 Gpc (right panels), 
corresponding to galactic halo and cosmological scenarios, respectively. The expansion opening 
half-angle is set to be Gb = 90° , and the variability timescale is At = 1 ms. Flat portions of 
the curves define the thin-shell limit AR/Rq ^1-/3, while the complementary sloping portions 
correspond to thick-shell expansions. 
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Fig. 4. — The minimum bulk Lorentz factor T^in for the most recent three of the six EGRET 
bursts depicted in Fig. but now for At = 1 second and d = 100 kpc (left hand panel) and 
d = 1 Gpc (right hand side), to illustrate the effect of lengthening At . Again, the curves represent 
T^^{et) = 1 solutions to Eq. (p7|), and source parameters are taken from Table |2[ Comparison with 
the bottom panels of Fig. ^ reveals that the Fmin are reduced significantly for these much longer 
variability timescales, however the dependence on At is nevertheless somewhat weak. 
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Fig. 5. — The minimum bulk Lorentz factor Fmin for the most recent three of the six EGRET 
bursts depicted in Fig. ^, for At = 1 second, but now as a function of the expansion opening 
angle 0b • The left hand panel is for d = 100 kpc and AR/Rq = and the right hand 

side depicts the cosmological case of d = 1 Gpc and AR/Rq = 10~^ ; the choices of AR/Rq 
correspond to transitions between the thin and thick-shell regimes. Again, the curves represent 
Tyy{et) = 1 solutions to Eq. (|37|), and source parameters are taken from Table ^ Clearly evident 
is the independence of F^jn on opening angle when 0^ » AR/Rq , and a very weak power-law 
dependence [whose slope depends on a burst's spectral index a via Eq. (|4^)] for ©b ^ AR/Rq . 



-49- 




Fig. 6. — The phase space diagram for the observational parameters Gmax and a , consisting of 
contours of constant T , as labehed. The contours correspond to solutions of r-^-y = 1 for specific 
choices of F ( = Fmin ) typical of both galactic halo sources at d = 100 kpc (three solid lines for 
each of F = 10,20) and cosmological scenarios at d = 1 Gpc (dashed lines: F = 300,600). The 
heavy contours represent fractional thicknesses 6 = AR/Rq = 1-/3, and for the galactic halo cases 
of F = 10, 20 , alternative thicknesses are represented by the light contours, with the lower curves 
for each F denoting 5 = 0.3(1-/3) and the upper ones 5 = 3(1-/3) . For all contours, the opening 
angle was set at 0b = 90° , the variability timescale set at 1 ms, and the flux /(I MeV) at 1 MeV 
assumed a value typical of EGRET sources. The observed EGRET values of ^max and a for six 
bursts (see Table §) are plotted as points; these bursts all have values of /(I MeV) different from 
the "canonical value" chosen here. 
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Fig. 7. — The phase space diagram for the theoretical parameters 0b and b = AR/Rq , consisting 
of contours of constant T , as labelled. As in Fig. ^, the contours denote solutions of r^-y = 1 for 
specific choices of F typical of both galactic halo sources at d = 100 kpc (three solid lines with 
F = 10,20,40) and cosmological scenarios ai d = 1 Gpc (dashed lines: F = 250,500,1000). For 
all contours, the maximum energy and spectral index were set at Smax = 100 MeV and a = 2.21 , 
respectively, and the variability timescale At and the flux /(I MeV) at 1 MeV are representative 
of EGRET-like sources. The vertical portions of the F = 10 and F = 250 contours define thin- 
shell/thick-shell regimes where the optical depth is independent of 0b • The three light dotted 
lines are the asymptotic limits given in Eq. (E^). 
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Fig. 8. — Solutions to Eq. (47) for the minimum bulk Lorentz factor T^in that guarantees source 
transparency up to energy Emax for two different source distances d as labelled. Values of Smax 
are chosen to probe beyond the EGRET energy range up to the domain of air Cerenkov detection 
techniques. Again infinite power-law source spectra are assumed. The source flux / at 511 keV 
is typical of BATSE burst detections; for this flux, the MILAGRO experiment will be sensitive to 
bursts with ah ^ 2.6 . The cosmological redshift was taken to be z = for simplicity. 



